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ABSTRACT 


The  primary  objective  of  this  study  was  to  provide 
a  contribution  to  the  phenomenon  "Suppression"  as  an 
aspect  within  the  military  environment. 

Analytical  models  explaining  these  aspects  were 
developed  in  order  to  identify  the  influences  to 
suppression.  Techniques  are  examined  for  including  the 
suppressive  effects  of  weapon  systems  in  Lanchester 
type  combat  models,  which  may  be  useful  in  wargame 
evaluations  of  military  judgements,  and  in  force  level 
planning.  The  study  also  provides  techniques  to 
analyze  and  fit  experimental  data  to  the  analytical 
models. 

The  data  to  verify  the  models  were  obtained  from 
related  experiments  performed  by  Combat  Development 
Experimentation   Command  (CDEC)  ,  Fort  Ord,  California. 

The  result  for  the  modelling  approach  to 
suppression  indicates  source  dependences  on 
guantitative  as  well  as  on  qualitative  features. 

The  functions  are  left  quite  general,  although 
some  functional  forms  are  derived  and  discussed. 
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I.   FORMULATION  OF  THE  PROELEM 


A.    INTRODUCTION 


In  recent  years  considerations  of  Human  Factors  aspects 
in  the  military  has  gained  more  and  mora  importance.  The 
study  of  psychological,  physiological,  and  environmental 
conditions  and  their  influences  on  the  performance  of  men  in 
"Man-Machine-Systems"  supports  the  development  of  new 
doctrines,  design  of  weapon  systems  as  well  as  training 
programs  for  troops. 

One  phenomenon  in  military  man-machine-systems  is 
suppression.  Modelling  this  phenomenon  has  only  recently 
been  given  much  attention. 

This  is  partly  because  modellers  did  not  understand 
exactly  neither  the  causes  of  suppression  nor  how  it  affects 
the  course  of  combats. 

There  is  an  intuitive  feeling  that  when  a  soldier  or  a 
combat  unit  is  being  fired  upon,  it  will  be  less  effective 
than  when  it  is  not  receiving  fire.  This  is  generally 
referred  to  as  suppression,  but  it  can  include  much  more.  In 
the  broadest  sense,  suppression  can  effect  individuals, 
units,  or  weapon  systems  in  different  types  cf  combat. 

This  paper  will  limit  itself  to  individuals  or  small 
infantry  units. 


In  order  to  be  able  to  build  up  functional  relationships 
based  on  the  suppression  idea_,  between  different  categories 
of  people  and  time,  this  paper  will  first  examine  different 
rates  as  input  to  analytical  Lanchester-type-models .  The 
rates  will  be  based  on  stochastic  aspects.  The  modelling  is 
done  under  different  viewpoints;  their  results  will  be 
compared  and  discussed. 

In   the  second  part  of  the  thesis,  the  parameters  cf  the 

models  will  te  tested  in  a  regression  analysis  against   real 

world   data  which   were   placed  at  the  writers'  disposal  ay 
CDEC. 

The  goal  of  this  paper  is  not  tc  present  a  final 
framework  cf  suppression;  however  it  may  contribute  to 
clarification  of  aspects  of  suppression,  and  help  to  embed 
it  as  a  component  in  future  large  size  models. 


B. 


BACKGROUND 


The  strains  imposed  on  individuals  in  our  society  are 
constantly  increasing.  Modern  technologies  and  constant 
efforts  for  improvements  of  standards  of  living  lead  to 
growing  difficulties  in  adjustments  or  sometimes  to  complete 
failures  to  adjust.  These  facts  create  stress  and  ve  may 
observe  that  the  degree  of  stress  increases  with  the 
difficulty  of  the  adjustment-problem. 


The  term  stress  will  be  used  as  a  substitute  for  what 
might  be  called  otherwise  as  anxiety,  conflict,  emotional 
distress,  extreme  environmental  conditions,  ego-threat, 
frustration,  threat  to  security,  tension  or  probably 
arousal.  Stress  can  be  thought  of  as  the  result  of  almost 
any  environmental  interference  [ Applay , p. 2  ]. 


The  stress-generating  features  of  the  civilian 
environment  are  great,  but  the  environment  created  by  modern 
warfare  possesses  additional  features  which  result  in  an 
increase  of  stress.  The  combat  environment  created  fcy  the 
weapon  power  of  the  enemy  causes  a  constant  threat  to  life. 
The  soldier  has  to  operate  under  this  threat  and  naturally 
he  will  respond  with  constantly  recurring  fear.  This  may 
break  down  the  soldier's  psychological  and  physiological 
resistance. 

Fear  and  anxiety  in  battle  is  common,  being  experienced 
by  between  80  and  90  %  of  combatants.  Pains  in  the  stomach, 
fatigue,  dizziness,  perspiration,  and  enhanced  heart-beat 
are  some  vegetative  correlates  to  fear  and  anxiety.  Of 
course  the  moment  when  the  individual  soldier  reaches  his 
breaking  point  varies  and  depends  en  individual 
psychological  and  physiological  resistance  and  the  severity 
of  the  battle. 

One  interesting  observation  from  Agiell  was  that 
auditory  sensations  convey  the  stress  of  battle  most 
strongly  and  most  directly.  The  psychological  effect  of 
weapons  is  directly  related  tc  their  sound  level  and  the 
freguency  with  which  the  sound  occurs  [ Agrell ,p. 2 1 5  ]. 

Each  enemy  grenade  causes  the  soldier  to  react 
constantly  with  fear.  Stress  and  fear  can  have  a 
significant  sensory  operating  characteristic,  e.  g.  the 
detection  threshold  and/or  sensitivity  may  decrease  as  a 
result  of  stress  [Weltmann,  G,  Christianson,  R.A.  and 
Egstrom,  G.H.,  p. 423-430]. 


C.   MILITARY.  SUPPRESSION  IN  COMBAT-ENVIRON MENI 


When  weapons  are  used  in  combat,  there  are  two  types  of 
effects  that  they  have.  The  first  type  of  effect  is  physical 
damage  or  injury  to  the  target  and  the  second  type  of  effect 
is  psychological. 

This  second  effect  of  a  weapon  has  led  to  the  term  or 
concept  "Coabat  Suppression".  There  exist  a  general  belief 
that  fire  suppression  is  important,  but  the  importance  of 
suppression  effect  on  combat  outcomes  as  compared  to  the 
effect  cf  other  areas  such  as  firepower,  mobility, 
intelligence,  command/control  has  not  bean  quantified 
adequately. 

1.   Definitions 


Suppression  can  be  generally  defined  as 
the  temporary  degradation  in  the  quality 
of  performance  of  an  individual  soldier 
or  unit  by  an  internal  or  external 
stimulus. 

CDEC,  1977 

A   more   useful   operational  definition   in   terms   of 

performance  capability  changes  is  provided  by  the  "Report  of 

the   Army   Scientific   Adviscry  panel  Ad  Hoc  Group  on  Fire 
Suppression" . 
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Ad  Hoc  Group,  1976 
This   definition   emphasizes   that  suppcession  is  not  a 
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single  effect  which  can  be  measured  totally  on  a  single 
quantitive  scale.  Suppression  effects  are  multidimensional 
and  the  "amount  of  suppression"  varies  among  these 
dimensions  (e.g.  fire  impact  points,  soldier's 
characterists,  and  reaction  to  the  fire,  his  combat 
experience  etc. )  . 

Suppressive  fire  in  a  combat  environment  can  suppress  a 
number  of  combat  activities;  for  example:  firing,  search  for 
and  observation  of  targets,  movements  of  units  or  command 
and  control. 

2-   Structure  Of  Fire  Suj:j;ression  Process 

Fire  suppression  is  a  complicated  process  involving 
many  physical,  environmental,  physiological,  behavioral,  and 
operational  variables.  The  important  point  to  emphasize  is 
that  the  behavior  involved  is  in  response  to  stimuli  that 
originate  both  externally  (combat  environment)  and 
internally  (personal  background,  training  and  experience)  to 
the  soldier  suppressee.  These  aspects,  however,  are  not 
included  in  this  paper. 

The  intensity  and  duration  of  suppression  can  not  be 
predicted  from  a  knowledge  of  the  combat  environment  alone. 
It  requires  an  analysis  of  the  underlying  motivational  and 
cultural  factors  and  of  the  context  of  the  combat 
environment . 

Xha  Fig  1  [Ad  Hoc  Group,  p.  36 ]  shous  a  schematic 
description  of  a  process  when  suppressive  fire  is  delivered 
and  its  affect  on  the  combat  result. 
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SCHEMATIC    STRUCTURE    OF    THE    FIRE    SUPPRESSION 
PROCESS      (Ad  Hoc  Group,    p.   56) 
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The  process  is  discribed  by  individual  functions. 


a.   Signal  Process 

Ihe  first  process  in  fire  suppression  is  the  generation  of 
signals  provided  by  suppressing  weapons. 

Inputs  to  this  process  are  the 

characteristics  of  the  weapon  systems  (caliber, 
amount  of  porpellant,  warheadtypes  etc)  and 

the  environmental  characteristics  (trajectory, 
platform,  arrival  points  etc). 

Characteristics  of   the  weapon  systems  can  vary  in  order  to 
increase  or  decrease  suppressive  effects. 

Some   parameters   that   are   considered   to   be  important  to 
suppression  signals  are: 

Muzzle  velocity  (an  increase  in  muzzle  velocity  is 
associated  with  an  increase  in  signal  variables)  ; 

Caliber  (as  caliber  increases,  the  firing  signals 
and  projectile  signals  increase  along  with 
lethality) ; 

Projectile  weight  (  penetration  depends  on  weight 
and  velocity  at  impact  and  increases  shock  coupling 
to  ground)  ; 

Warhead  charge  weight  (  the  explosive  charge  weight 
determines  the  energy  in  the  pressure  pulse)  ; 

Additional  parameters  like  fire  freguency  and 
proximity  of  shots  could  also  be  mentioned. 

Environmental   characteristics   are   also  variaoles  but  thay 
can  not  be  determined  completely. 

Environment   has   an   influence   on   the   signal 
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generation  -and  transmission-process.  For  example,  auditory 
signals  that  result  from  the  impact  of  projectiles  depend 
heavily  on  the  nature  of  the  object  or  material  impacted. 

A  soft  yielding  material  such  as  dusty  ground  or 
sand  receiving  the  impact  of  a  projectile  will  produce  a 
different  pulse  and  sound  than  will  hard  unyielding  ground 
under  the  sane  impact. 

Sound  signals  can  be  attenuated  by  the  shadowing 
effect  of  large  obstacles  or  may  be  increased  by  echo  or 
reverberation.  Visual  signals  are  strongly  modified  ty  the 
condition  of  lighting.  Haze,  fog,  rain,  and  snow  act 
similarly  to  smoke  and  dust.  The  visual  field  is  also 
reduced  and  interrupted  by  terrain  and  other  obstacles. 

The  environment  modifies  the  produced  stimuli 
when  they  are  transmitted  to  the  location  of  the  suppressee. 
Outputs  of  this  process  are  the  attenuated  sensory  signals 
that  become  imputs  to  the  hunan  sensory  receptors. 

b.   Human  Process 

Many  of  the  determinants  of  the  soldier's  performance  on 
the  battle  field  are  unkown  or  at  least  uncertain  -  thought 
of  as  influenced  oy  chance  factors.  This  emphasizes  the 
difficulties  of  predicting  human  benavior  in  a  combat 
environment.  Tne  human  process  (sensory  and  perception) 
converts  the  received  signals  into  a  perception  of  the  risk. 

Eattle field  stimuli  effecting  the  individual  are 
detected  and  converted  into  sensory  data  by  a  process  such 
as  vision  and  audition,  so  the  sensory  process  suggests  that 
the  weapon  systems  stimuli  relevant  to  suppression  are  the 

-loudness,  and 

-visual  impact. 

There  exist  moderating  factors  that  influence 
the  operating  cnaracteristcs  of  the  sensory  process  and  that 
determines  which  stimuli  are  effective. 
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Sensory  modifiers  (i.e.  earplugs, 
night-vision-devices)  serve  to  change  users'  sensitivity 
range.  A  major  effect  of  these  devices  is  to  change  the 
salience  of  stimuli. 

High  concentration  on  an  activity  or  a  high 
level  of  effort  on  an  activity  (e.g.  missile-gunner  is 
tracking  a  target  or  reloading  his  system)  may  increa-se  the 
absolute  threshold. 

Ihe  posture  of  a  soldier  (standing  or  sitting) 
and  the  sequence  of  posture  (observing,  ducking,  observing) 
influences  the  sensory  capabilities  (e.g.  observing  for  10 
seconds  continuously  is  not  equivalent  to  observing  5 
seconds  ,  ducking  10  seconds  and  observing  5  seconds). 

Ihe  perception  process  integrates  sensory  and 
other  information  into  a  perception  of  the  risk.  Risk  refers 
tc  the  uncertainty  of  damage,  injury,  or  loss.  It 
characterizes  decision  situations  in  which  the  consequences 
of  choosing  an  action  are  uncertain. 

If  there  is  no  uncertainty  in  the  possible 
outcome,  there  is  no  risk. 

Perceived  risk  is  a  function  of  uncertainty  and 
the  subjective  value  the  individual  associates  with  each 
outcome. 

Eerceived  risk  represents  the  output  of  the 
combined  'sensory  and  perception  process.  It  depends  en  the 
individual's  experience  and  training  in  assessing  risk  from 
sensory  information.  Also  cover  provided  by  the  environment 
and  the  individual's  pesture  may  influence  risk   perception. 

c.   Reaction  Process 

Given  the  input  perception  of   risk,   this   process   causes 
physical  and  mental  reactions,  which  depend  on  the 
-current  mission 
-task 
-activity 
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-combat  training  doctrine  and  experience 

-group  dynamics  and 

-the  quality  of  leadership. 

It  is  conjectured  that  two  individuals  who  perceive  the  same 
degree  of  high  risk  ,  but  who  have  different  amounts  of 
ccmbat  engagement  experience,  might  be  lively  to  react 
differently  tc  the  risks. 

The  soldier's  reaction  is  also  influenced  hy  his 
current  state.  A  soldier  who  has  recently  ducked  may  be  more 
likely  to  duck  than  one  who  has  not  ,  given  the  same 
delivered  fire. 

Prior  reaction  or  sequence  of  reactions  may  te  a 
good  predictor  of  the  coming  reaction. 

d.   Performance  Effects  Process 

Given  the  reactions  of  the  human  behavior  process,  it  is 
conjectured  that  these  directly  affect,  the  performance  of 
certain  activities  of  the  suppressee  in  a  calculable  way. 

If  for  example  the  suppressee  takes  cover,  he 
may  fire  less  often  and  less  accurately  and  also  might  be 
less  vulnerable.  The  magnitude  and  duration  of  these  changes 
in  performance  are  dependent  on  the  characteristics  of  the 
system  employed  by  the  suppressee  and  the  target  of  his 
activity . 

So  the  nature  and  duration  cf  change  in 
performance  capabilities  is  determined  by  the  performance 
effects  process. 

3  .   S  ujD  _p_ression  in.  Field  Experimentations 

A  fire  suppression  research  program  reouires 
significant  experimentation  on  behavioral  attitudes  and 
reactions  tc  risk.  This  necessity  causes  tremendous 
difficulties  in  trying  to  induce  actual  behavior  in  soldiers 
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in  field  experiments.  Former  studies  shewed,  that  the 
soldiers  felt  true  psychological  stress  only  in  situations 
in  which  they  believed  that  th-ey  were  in  real  danger.  Such 
situations  are  difficult  to  contrive  and  tc  control.  Social 
and  ethical  limits  and  legal  regulations  preclude  the 
introduction  of  actual  physical  risk.  Scldiers  must  be 
taught  the  "rules  and  risk"  defined  in  that  context'.  The 
success  of  playing  the  role,  being  an  individual 
participating  in  a  combat  engagement,  depends  on  the 
soldiers'  motivation  and  willingness. 

Because  of  these  reasons  and  the  nultidimensional 
shape  of  the  fire  suppression  process  as  mentioned  earlier, 
suppression  in  field  experimentations  may  be  restricted  only 
to  some  variables  involved  in  this  process. 


D.    APPROPRIATE  OBJECTIVES 

The  overall  objective  of  a  fire  suppression  research 
should  be  to  relate  changes  in  performance  capabilities 
caused  by  fire  suppression  [Ad  Hoc  Group,  p.  110]. 
Responding  more  directly  the  following  objectives  may  be 
determined: 

Indicating  the  effects  of  suppression  on  combat 
results,  i.  e.  to  develop  rates  cf  suppression. 
These  values  may  be  compared  to  other  effected  areas 
and  probably  employed  in  computer  sinulations.  These 
numbers  represent  two  kinds  of  variables:  Weapon 
system  variables  and  human  suppression  performance, 
given  operational  and  environmental  conditions. 

Determining  characteristics  cf  suppressive  fire 
systems,  characteristics  which  should  be  assigned  to 
such  a  weapon  system.  Results  are  developed 
experimentally  .  Chapter  III  B.  will  display  some 
evaluated  parameters  and  constants  for  the  developed 
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model,  which  may  also  represent  the  suppressive 
characteristics  of  the  used  weapon  systems. 

Heducing  suppressive  effects.  Ways  to  reduce  the 
effect  of  suppressive  fire  may  also  be  considered  as 
an  appropriate  objective  of  suppression  research. 
Special  training  or  equipment  can  be  assigned  to  tne 
soldiers  or  new  tactics  can  be  developed.  This 
objective  is  beyond  of  the  research  of   this   paper. 


E. 


POSSIBLi  ALTERNATIVES 


In  ordei  to  get  information  anouc  the  fire  suppression 
process,  previous  investigations  were  based  on  interviews 
and  questionaires ,  because  valuable  information  of  the  fire 
suppression  process  is  stored  in  the  minds  of  combat 
veterans.  Studies  on  veterans  of  the  Vietnam  conflict  and 
the  wars  cf  the  Near  East  would  be  especially  useful,  since 
newer  weapons  were  employed  and  the  combats  were  shorter  and 
mere  intensive. 

These  studies  may  provide  a  good  insight  tc  the 
suppression  process  and/or  may  also  deliver  valuable  inputs 
for  the  modelling  approach. 
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II.   CONSTRUCTION  OF  AN  ANALYTICAL  MODEL 


A. 


ASSUMPTIONS  AND  SIMPLIFICATIONS 


Based  on  the  foregoing  discussion,  the  following  model 
is  a  detailed  model  [ Taylor , 1978, p .  12  ]  which  starts  cut  by- 
considering  the  behavior  aspect  of  a  human  being  under  the 
influence  of  artillery  fire  power.  It  is  assumed  that  such 
simple  uicdels,  which  represent  a  small  part  of  a  total 
scenario,  can  be  used  profitably  to  investigate  system 
dynamics  of  lore  complex  models.  The  value  of  the  model 
derives  from  the  fact  that  it  forms  intuitively  plausible 
and  transparent  subsets  in  a  large  composition  of  other 
subsets  which  determine  the  basic  structure  of  the  complex 
operational  model.  In  other  words,  the  whole  is  described 
in  terms  of  the  sum  of  its  parts. 

The  basic  concern  of  the  analytical  model  developed  here 
will  be  to  model  the  behavior  of  an  individual  experiencing 
artillery  fire,  considered  as  a  function  of  time,  where  that 
behavior  depends  upon  ammunition  types  and  location  of 
detonating  rounds. 

On  the  basis  of  particular  assumptions  and 
simplifications  it  will  be  possible  to  apply  the  results 
cbtained  to  a  group  of  people  (a  force)  en  a  battelfield. 
The  result  of  these  considerations  will  provide  a 
relationship  between  time  and  tne  actual  number  of  people 
affected  by  the  fire  power  of  the  artillery  .  This  last  step 
of  the  model  is  carried  out  by  using  Lanchester  type 
equations,  so  called  after  the  pioneering  work  of  F.  W. 
Lanchester.  Finally,  the  models  enable  one  to  estimate  the 
total  firepower  of  the  force  at  any  point  in  time;  depletion 
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of   total    firepower   is   caused   by   attrition   and   by 
suppression. 


E.   RATES  FOE  THE  MODEL 


1  •   Rates  Of  Suppression 


The  tasic  considerations  in  the  preceeding  paragraph 
support  the  assumptions,  that  the  behavior  of  a  suppressee 
can  be  expressed  by  a  conditional  probability  of  suppression 
as  a  function  of  miss  distance  r  and  aspect  angle  9  .  This 
function  is  represented  by  the  family  of  surfaces  of  the 
form: 


P(S/6,r)  =  exp[-  ~   r2(l  -  e  cos  9)]  (2.1) 


where  P(S/9,r)  is  the  conditional  probability  of  suppression 
given  that  a  particular  round  impacts  under  a  certain  aspect 
angle  9  and  a  certain  miss  distance  r  away  from  the  foxhole. 

The  line  along  the  angle  9=0° is  identical  with  the 
line  of  sight.  It  is  the  main  direction  of  ccservation. 

lhe  constant  K  and  e  are  parameters,  which  are 
determined  by  the  experiment  itself  and  by  the  environ  mental 
conditions. 

Ihey  can  be  influenced  by  factors  as  discussed  in 
chapter  I  which  may  be  recalled  here  briefly. 

-type  of  ammunition 

-freguency  of  arrival  of  rounds 

-perceptual  damage 

-total  time  spent  in  the  foxnole  (learning  process) 
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-noise  appearance  of  the  rounds 

-flash  light  intensity  of  the  rounds 

-performance 

-personal  factors  like  age,  personal  condition,  etc. 

-degree  of  stress 

-motivation 

The  computational  evaluation  of  the  constant  K  and  z  will  be 
performed  in  chapter  III  B.  and  C.  In  particular  it  will  be 
important  to  determine  K  as  it  varies  with  different  types 
cf  ammunition. 

The  mathematical  conditions  for  the  two  parameters  K  and  e 
are : 

0  <  e    <  1 

K  >  0 

If  6  is  held  fixed  ,0°<8  <360°,  and  ?(S/9,r)  varies 
between  0<P  (S/$,r  )  <  1  we  obtain  a  family  of  functions,  which 
is  two-dimensional  and  shows  an  exponential  relationship 
between    P(S/9,r)    and    r. 

This  is  illustrated  in  the  following  figure,  where 
the  angle  9  is  held  fixed  at  C'and  180*.  The  parameters  X  and 
e  are    assumed    as    being    1500    and    0.7    respectively. 
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r/m  =  distance/meter 


P(S/e0  ,r) 


r/m 

for  G  «1S0 


foxholo 


I 1 1 1 1 1 1  B»- 

100  , 

r/m 

for  900° 


Figure    2   - 


1     2 
FUNCTION      P(S/eQ,r)    =   cxp[-  -  r    (1   -  c    cos   9)] 
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If  en  the  other  hand  P(S/'0,r)  is  fixed, 
0<P  (S/ 0/ r) < 1 ,  and  the  equation  is  solved  for  r,  we  will 
obtain  an  "egg-s haped"  function  with  iso-levels  of 
probability  of  suppression  P(S/8,r).  The  foxhole  is  located 
in  the  middle  of  the  coordinate  system.  Along  9=0°,  the  range 
r  is  a  max  for  a  certain  fixed  probability  cf  suppression  , 
while  at  8=180°,  r  is  a  min  for  the  same  probability. 

The  function  takes  the  form: 


-K  £n  P(S/0 ,r) 
1  -  e  cos  8 


(2.2) 


Ihe    following    graph    shows    the      family      of      functions      for      4 
representative   selective  probabilities    of   suppressions. 


P  (S/8/r)=0.1«i 


i  =  1,2,3,4 


(2.3) 


The   parameters  K  and  e   are  again  assumed  to  be  1500  and  0.7 
respectively. 
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r/ra     -     distanco/mater 


j u 


0  r/m 

e=270c 


r/n 

"   .Cor  0-1 80e 


Figure    3    - 


ISOLINEb       r      =   —c    Cos   0       '      i   =   1'2.3.4' 
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These  "egg-shaped"-f unctions  of  iso-probabilities  of 
suppression  simulate  the  reaction  of  human  beings  looking 
along  the  main  axis  6=0* in  a  very  simple  way.  The  functions 
take  into  account  the  visual  and  accoustical  perception 
resulting  from  any  given  detonation  of  a  round,  where  noise 
and  light  are  the  major  stimuli. 

One  could  think  of  other  functions  wich  simulate  tae 
behavior  of  an  antitank  gunner  exposed  tc  artillery  rounds 
for  instance: 


1  2 
P(S/0,r)  =  L-exp[-  -  r  (1  -  e  cos  8)]  (2.4) 


Ihis  type  of  function  allows  the  probability  of  suppression 
E(S/0,r)  to  te  smaller  than  one  at  its  maximal  value  and  has 
the  same  general  behavior  as  the  function  before.  As  a 
third  modification: 


P(S/0,r)    =  exP[-  ||r|(l  -  e    cos   6)]  (2.5) 


these  functions  have  the  disadvantage  that  they  have  a 
discontinuity  at  r=0.  The  integration  which  is  necessary  in 
the  following  derivation  is  more  difficult  than  the  chosen 
one. 

However,  the  crosscut  sections  of  these  functions 
are  not  "egg-shaped"  iso-f unctions  but  rather  simple  conic 
sections  (ellipses)  where  the  foxhole  is  located  in  the 
center  of  one  focus  point. 

The  paper  will  continue  with  the  function  first 
described  in  (2. 1) ,  because  of  its  simplicity  and  variety 
cf  applicaticn. 

In  order  to  derive  a  rate  of  suppression  for  the 
model,   it   is   necessary  to  evaluate  now,  in  a  second  step, 
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the  unconditional  probability  of  suppression  P (S)  for  any 
incoming  round,  no  matter  where  it  will  impact  around  the 
foxhole.  This  will  be  performed  by  matching  the  conditional 
probability  P(S/9,r)  with  the  area-hit-probatility  P(A). 

Consequently,  to  the  engagement  procedures  of  the 
artillery,  the  targets  which  shall  be  suppressed  by  the 
artillery  are  categorized  as  small  personal  targets.'  This 
implies  tnat  the  mean  point  of  impact  (tfPI)  of  a  given  set 
of  rounds  lies  on  the  target,  which  means  also,  that  the 
density  has  its  max  value  at  that  point.  According  to  the 
U.S.  Army  Field  Manual  FM  6-161-1  page  2-2  it  can  be 
ascertained  that  the  MPI-error  is  destributed  normally  with 
its  mean  at  the  aim-point.  Thus  it  will  be  hypothesised  that 
in  sucn  a  case  the  distribution  of  incoming  rounds  is 
bivariate  normal  with  parameters 


v2   =  0 


o .  =  a„  =  a 


(2.6) 


In  order  to  simplify  the  model,  it  is  proposed  that  the 
dispersion  of  rounds  expressed  in  the  standard  deviation  is 
equal  in  all  directions.   With 


p  =  0 


(2.7) 


no  correlation  is  assumed  between  horizontal  and  vertical 
deviation.  The  normality  is  also  preserved  if  more  than  one 
artillery  gun  is  shooting.  The  essential  change  which  has 
to  be  made  when  a  whole  artillery  unit  will  deliver  the 
rounds  will  be  the  value  of  the  standard  deviation  a 
In  addition,  a     is  determined  by  factors  like: 

dispersion 

the  fact  that  the  location  of  the  target  is  only 
estimated  and  an  artillery  unit  delivers  rounds  in  a 
fire  area. 
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type  of  ammunition 


conditions    of   the    weapon   systems 
ballistic  properties  of  the  rounds 


(precision) 


distance  to  target 

caliber 

wind    and   other    weather   conditions. 

It    follows    frcm    the    hivariate   normal   density      for      artillery 
hits: 

1 


f(x   ,x2)    = 


2      2 
2ir(l-p   )<J 


(2.8) 


(  1 

x  exp  - 


(  2d-p2)  L\ai 


Vri 


Xl   rl\ /X2~r2 


+ 


x2-r2 


2-i 


and    by    inserting      the      previous      mentioned      assumptions      and 
changing    to    pclar    coordinates: 


2  2  2 

r     =  x.   +  x2 


cos   9   =  — 

r 


(2.9) 


that  the  unconditional  probability  of  suppression   P(S)   can 
be  written  in  the  form: 


G=2-rr   r=°° 
P(S)  =  /     /   P(S/0,r)«f  (r,0)  do  r dr 
e=o   r=0 


(2.10) 


Eecall  that  an  area  element   in   polar   coordinates   can   be 
ex  pressed: 


r-dG-dr  =   d(-j-)    d0 


Ihe   conditional    probaoility   of   suppression   is 

1      2 
P(S/8,r)    =  expt-  -  r    (1   -   cos   9)] 
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(2.11) 


and   the    density   for   artillery   hits   f(r,9)     is 


f(r,8)  = 


2tto 


2      exPt"  I    <f>     ] 


(2.12) 


inserting  both  into  the  expression  anove  ,   we   can   perform 
the  integration: 


Q=2t7  j-=oo  j  2 

P(S)    =     /  /      exp[-  £  r2(l   -   c    cos   9)   — ^y  exp[-  j(~)    ]-d(y-)    d9 

6=0     r=0  2,a  (    } 


after  a  change  of  integration  variable; 

2 


2      n   =  0 
r=0 


dOy-)    =   dz 


=    z 


Z  =    oo 


(2.14) 


we    find 

PCS) 
p(s) 


9=2ir    z=c° 
=      /  J      exp 

6=0      z=0 


9=2tt 
/ 


2    ,-,  1 

-    (1   -   £    cos   6)    +  — 

o 


2ua' 


dz   d9 


0    ~2':'jZ  [\   0-  -  e    cos   6)    +  ~ 

a 


•  exp 


/2 


-  (  -   (1   -  £    cos   a)   +  — 


z=CJ 


de 


PCs)  = 


9=2tt 
j 1 

6=0      2tto2/  |-  Y  cos   8   +^2 

2Ka 


d0 


(2.15) 


K 


0=2t7 


PCS)  =-±j    / 


4tto      6=0       (1   +  -~)  +   C-e)    cos 
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•  de 


(2.16) 


by   replacing 


b  =  1  + 


_K_ 

2a2 


c  =  -e 


(2.17) 


the  integral  reduces  to  a  known  form  which  can  be  solved. 


P(s)  .  -JL.  T  *- 

/      2  Jn     b  +  c*cos 


2  2 

where      b      >   c 


k     J7         de 


p(s)  -ri/.bT 


.,      -   ,,    ,.,   b   +   ccos   6         „      2 
2tto      9=0  2ira 


(b-c)    tan 


9    -i 


arctzan 


/X7  vV^ 


b  -c 


PCS)  =-^         2 


2*°   yx? 


==>  p(s) 


2a    v  b   -c 


(2.18) 


Inserting   back    the    expressions    for    b   and   c,    one   obtains: 


PCS)    = 


(2.19) 


2o    \\    1  + 


-   c 


2a' 


After  the  data  analysis  in  chapter  III  B.  and  C. 
where  the  parameters  K,  e  and  0  will  be  determined,  it 
will      be    possible   to    evaluate   the    probability   of    suppression 
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E(S)  according  to  the  preceeding  formula. 

This  probability  adequately  models  those  cases  where 
temporary  suppression  is  the  only  possible  response  tc  a 
given  detonation  of  a  round.  It  may  be  useful  to  expand  the 
model  by  introducing  a  second  paired  outcome  at  any  given 
impacting  round.   Until  now,  we  had  one  paired  outcome: 

The  individual  was  either  suppressed 

cr  not  suppressed. 

In  addition  to  these  we  consider  a  second  paired  outcome: 

The   individual   is   permanently  suppressed  by  being 
wounded  or  killed 


or  net  permanently  suppressed. 

If  we  focus  the  attention  on  modelling  the  effects  under 
this  expansion,  it  is  possible  to  evaluate  a  new  probalility 
of  suppression  ?(S)  and  a  probability  of  kill  P  (K)  by  Baking 
the  following  assumptions: 

The  conditional  probability  of  kill  unlike  the 
ccndtional  probability  of  suppression  does  not 
depend  on  the  aspect  angle  9. 

It  can  be  represented  by  a  smooth  curve  of  the 
following  form: 


P(K/r)  =  exp(-  ~   r2) 


(2.20) 


where  H  is  a  positive  constant  i.e.  H>  0.  This  function  was 
selected  on  intuitive  grounds  rather  then  based  on 
real-world  data.  Taking  a  vertical  cut  through  the  surface 
function  above  along  any  angle  9,  one  obtains  the  following 
figure.  The  parameter  H  is  arbitrarily  chosen  as  being 
H=50.  Clearly  this  parameter  is  a  function  of  different 
factors  like  terrain  and  ammunition. 
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r/m  =  distance/meter 


P(K/r) 


r/m 

for  9-180° 


foxhole 


for  0O°O£ 


1     2. 


Figure -4  -     FUNCTION       P(K/r)   =  exp  (-  -  r   ) 
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If  on  the  other  hand  P  (K/r)  is  fixed,  0<P|K/r)<1,  and  the 
equation  is  solved  for  r,  one  will  obtain  concentric  circles 
with  iso-levels  of  probability  of  kill  P  (K/r)  around  the 
foxhole. 

The  function  has  the  form: 


r   =  -  H-2,n  P(K/r) 


<  (2.21) 


The    following      Fig    5    shows    the    family      of      functions       for      4 
representative   selective  probabilities   of    sup fressions. 


P  (K/r)  =0. 1  -i 


i=1,2,3,4 


(2.22) 


The  parameter  H  is  again  chosen  to  be  50. 


32 


r/m 
0-270° 


foxhole 


r/m 

for  0-180° 


r/m  ■  distance/meter 


Figure    5    -       ISOLINES      r     =  -H-£n  P(K/r) 
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The  concentricity  around  the  foxhole  expresses  the 
fact  that  no  matter  in  which  direction  the  person  in  the 
hole  is  looking  the  kill  effect  is  just  determined  ty  the 
distance.  It  would  be  beyond  the  paper  to  verify  this 
assumption  ,  and  it  would  be  extremely  difficult  to  collect 
data  for  it.  The  reader  must  be  content  for  the  moment  with 
the  earlier  presented  intuitive  argument. 

This  newly  introduced  function  leads  to  a  revision 
cf  the  probability  of  suppression  P(S)  and  the  evaluation  of 
the  unconditional  probability  of  kill,  P(K).  The  density 
assumed  earlier  for  artillery  hits  is  used  again: 


f(9,r)  =  j  exP 

2ira 


2  V 


(2.12) 


Probability    cf   suppression    P  (S)     but   not    kill    is: 

9=2tt    r=°° 
P„(S)    =     /  /      P(S/0,r)[l   -   P(K/r)]-f (0,r)    r   dr   dG      (2.23) 

6=0      r=0 


Remark:  The  subscript  K  at  p  is  used  to  indicate  that  this 
suppression  probability  appears  together  with  the 
probability  cf  kill  P(K). 


The   expression   ?  (S/9  ,r).(1  -P  (K/r)  )   means   that    a 
suppressed  but  did  not  kill  the  individual  at  (r,9) ► 


round 


Erobability    of    kill    P(K) 


9=2tt   r->-°° 
P(K)    =  /  /        P(K/r)-f(e,r)    r  dr  d9 

9=0      r=0 


(2.24) 


Ihe  computation  for  E(S)  and  P(K)  runs  along  similar  lines 
as  in  the  earlier  evaluation  of  P(S).  Fci  the  specific 
functions  we  find: 


3a 


0=2-rr   r=<»  i      ?  l  i  2  2 

P(K)    =     /  /     exp(-  i  r   )    •   —4-  exp[-  j  (J)    ]    d(y-)    dO 

0=0      r=0  2tto 


(2.25) 


9=2tt   r=°°  ,      2  12 

P|/s)  =  /  /     exP["  7  r    (1  -  g    cos   0)]    [1   -   exp(-  -  r  )  ] 

K  6=0      r=0  *  H 


2ira 


1      r   2  r2 

exp[-  ±   £)     ]    d(-y)    dG 


(2.26) 


change  again  the  variable  of  integration 


=   z :  z 

'  r  -*■ 


£  •*-  oo 


zr=0   ■   0;  d(-y-)    =   dz 


(2.14) 


9=2      z  +  » 


P(K)    =     /  / 

P(K)    = 


exp 


0      z=0      2ttcj 
1 


a 


dz   d0 


(a    /H)    +  1) 


(2.27) 


(2.28) 


0=2tt     2-*»oo 


PK(S)    =     /  /      exp  -[£  (1   _  e    C0:3   9)    +    1  j2   |        1       dz   d 


0  (2.29) 


0=0      z=0 


iva 


"     /  /        T  exP  Mv    (1   -  c    cos   0)    +  -  +  "V]z'    dz   d6(2.30) 

0=0      z=0      2tto  '      K 


H  2 

a 


Recognizing      that      the      first   double    integral    is   exactly    the 
same   as    fcefore   we   can   simplify    to 


PK(S)  - 


2.2\/(l  +     K 


2a 


-    E 


)=2tt 

/  TfT 

3=0      2™      K 


cos    0 


2         1 

+  -  +  -^ 
H  2 

a 


d6 


(2.31) 
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PR(S)    =  P(S)   - 


K 


0=2tt 

7         /  / 

4tto       6=0 


dO 


1   +  -\  +  |  )  +■    (-:)    cos 


2a 


(2.32) 


The  integration  can  be  done  by  using  the  same  formula  as  for 
the  computation  of  P(S)  ,  except  that  this  tine  b  is: 

b»(1  +  JL  +  in  (2.33) 


2a 


2   H 


From  this  it  follows: 


PK(S)  =  P(S)  - 


Vs)  =  7T 

2a 


2a 


1  + 


K 


2a' 


—  z 


(2.34) 


(2.35) 


After  having  derived  the  desired  probaoilities  as  summarized 
on  the  following  figure,  we  are  able  tc  evaluate  the 
different  rates  of  temporary  and  permanent  suppression. 

Ey  assuming  a  certain  fire  rate  A  with  which  the 
artillery  is  firing  in  the  area  where  tha  anti-tan/:  gunners 
are  located,  the  different  rates  will  have  the  following 
form. 


iwo-event  model  (unsuppressed-suppressed) 

Bate  of  suppression  X     =  X^   P(s) 
r  s    f 


(2.36) 
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Four- event- mod  el     (  un  sup  pressed- sup  pressed-  sue  vived-lcil  led) 

Rate    of    suppression      X     =  \   •?   (s)  (2.37) 

S  1  t\ 


Hate    cf   killing 


Xk  =  Xf  P(K) 


(2.38) 


One  important  rate  oust  still  be  developed;  it  is  the  rate 
cf  rise  from  a  suppressed  state.  He  model  this  by  stochastic 
process  metheds. 
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PCS)  - 


2o2  ,1*4 


2a' 


-  e 


(2.19) 


P(K)  = 


H 


+  1 


(2.28) 


PK(S)  -  P(S)  - 


-2V(i+^r-*2 


(2.34) 


Figure  6  -   SUMMARY  OP  THE  DERIVED  FORMULAS 
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2.   Bate  Of  Rise 

The  model  so  far  represents  a  situation  in  which 
there  are  three  possitle  outcomes.  When  any  independent 
round  impacts,  the  person  in  the  foxhole  is  assumed  to  be 
either 

suppressed, 

net  suppressed,  or 

is  killed. 

Naturally  this  is  only  true  if  we  assume  that  he  is  always 
hack  up  again  when  the  next  round  impacts.  This  fact  trings 
up  the  necessity  of  considering  the  process  along  the  time 
axis. 

Eecause  of  the  complexity  of  a  stochastic  process 
which  reflects  the  behavior  situation  of  an  individual 
exposed  to  arriving  artillery  rounds,  it  seems  useful  to 
start  with  a  very  simple  process,  in  which  the  time  of 
suppression  fcy  a  particular  round  is  considered  to  be  fixed. 
In  a  next  step,  this  constant  response  time  is  randomized 
over  the  numcer  of  people  under  consideration,  i.e.  each 
individual  has  his  own  random  time  which  however  is  assumed 
constant  for  the  process  itself. 

Finally,  the  response  time  of  a  single  person  say  be 
considered  to  be  a  random  variable  coming  from  a  certain 
distribution. 

Eoth  approaches  are  simple,  and  represent  a  first  atteupt  to 
describe  the  actual  situation  in  different  ways.  Certainly 
the  suppression  time  can  also  be  considered  as  a  function  of 
both  above  mentioned  processes  or  of  the  miss  distance  r  and 
the  aspect  angle  9,  or  of  a  combination  of  all.  But  these 
dependencies  would  rather  complicate  the  mathematical 
derivations  and  lead  beyond  this  study. 

Both   processes   consider   rounds  that  arrive  in  the 
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neighborhood  of  the  foxhole  in  accordance  with  a  Poisson 
Erocess  having  a  rate  ^f  .  The  model  consists  of  two  process 
states.  One  is  the  "up"-state  which  means  that  the  person  in 
the  hole  is  unsuppressed  and  able  to  act  according  to  his 
mission.  The  other  state  is  termed  the  "down"-state  which 
results  from  an  incomming  round  and  a  possible  reaction  to 
it  in  the  form  of  suppression.  In  this  state  the  perscn  is 
not  able  to  fullfill  his  mission,  he  is  physically  down. 

The  first  process  assumes  that  the  person  in  the  foxhole 
changes  froi  the  "down"-state  to  the  "up"-state  only  if  he 
recognizes  a  gap  of  at  least  T  time  units  before  the 
appearance  of  the  next  round.  This  means  that  the 
"down"-'State  has  a  duration  of  exactly  time  T  if  no  rounds 
arrive  in  the  time  interval  (0,T) .  The  time  T  in  connection 
with  such  a  process  may  fce  called  the  critical  gap. 
[Gaver, p . 48 1  ].  The  following  figure  shows  the  basic 
relationship : 
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FIRST  CASS  :   t'  >  T 


State  0  "up" 
State  1   "down" 

T  =  critical  gap 


A 


State 


Impact  of  first 

suppressing 

round 


Impact  of  second  suppressing 
round   (process  starts  all 


over  again) 


t/sec 


5C0ITD  CASS  :   t'  ^  T 


A  State 


0  0 


Impact  of  first 

suppressing 

round 

1 • 

T 


Impact  of  second  suppressing 
round  (process  starts  all 
over  again) 


1 


^■ 


-J 


t/sec 


Figure  7  -   SCHEMATIC  DESCRIPTION  OF  IHE  MODEL 
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For  the  time  being  T  is  assumed  fixed.   However  it  is   mere 

reasonable   to   describe   it  by  a  random  variable,  since  the 

time  gap   T   is   determined   by   factors    like  learning, 

accustoming  to,  cr  overcomming  of,  fear,  stress  etc.  The 
assumption  of  a  fixed  T  has  to  be  changed  in  a  later  step, 
where  it  will  be  defined  as  a  random  variable. 

Furthermore,  the  change  from  the  "up"-state   to   the 

"down"-state  is  accomplished  by  arriving  rounds  to  which  the 

person  reacts  through  suppression  ("down"-state)  with  tne 
earlier  computed  probability  of  suppression  P  (S)  . 

Now  let  t  be   equal   to   the   total   time  of   being 

continually  suppressed.  It  is  possible  tc  define  x  in  the 
following  way: 

!T  if  no  suppressing  rcund  arrives  in  (0,T) 

t'  +  x'  if  the  next  suppressing  round  falls  (2.39) 
before  T. 

In  this  case  the  random  variable  t'  is  an  independent  version 
cf  the  random  variable  t. 

According  to  this  definition  it  is  possible  to  compute  the 
expected  value  of  x  which  represents  the  mean  time  spent  in 
the  "down"-state,  or  the  mean  time  of  being  continually 
suppressed.  The  symbol  used  for  it  will  be  E[t]. 

In  order  to  get  the  arrival  rate  cf  suppressing 
rounds  which  contribute  to  the  time  of  being  suppressed,  we 
have  to  multiply  the  fire  rate  A-  with  the  probability  of 
suppression,  P (S)  ,  which  was  evaluated  in  chapter  II.B.1.  We 
cbtain 

As  =  Af-P(S)  (2.36) 

First    we    have    to   state    the    following    two    expectations: 


E[x|no   round   in     T]    =  T 


E[t I  round   In     T]    =   t' 


(2.40) 
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where  t' =  time  of  second  suppressing  round.  Removing  the 
condition  on  t'  leads  to  expected  time  of  being  suppressed 
under  the  influence  of  incomming  rounds. 


-\   T         T     -X    t'  -XT 

E[T]    -  T«e      S    +    J      e      S     X    dt'ff    +  E[t]]+B[t]    [1-e      S    ] 
0  S 


E[t]    = 


-X    T 
T.e     S     +  X 


-     -X    t 

s 

- 

e 

,  2 

(~ 

t 

"       1) 

X 

—                     <5 

_ 

T   T 


->0 


-X   T 
+  E[t]    (1-e      S   ) 


(2.41) 


-XT.  -XT  -XT 

E[t]    =Te      S     +  f  [1  -    (1  +  X   T)    e      S]+  E[t](1   -   e      S    ) 

A  S 


.-XT 
E[t]    =  T    +         _x   T    [1  -    (1  +  XsT)e      S    ] 


E[t]    =  r±-   (e   S      -  1) 


(2.42) 


inserting     X     =  X   'P(S)       we  obtain 

s    r 


,      X  p(S)*T 

E[T]=i7(sy(e  -« 


(2.43) 


Extending  the  idea  for  a  fixed  time  gap  by  applying 
it  to  a  group  of  people  separately,  we  are  able  to 
reformulate  the  process  in  the  following  way: 


Assume  a  group  of  perscns  in  which  each  individual  sticks  to 
a  certain  but  randoi  time  gap  T  when  responding  to  a 
suppressive  round,  and  assume  that  this  random  variable  T 
based  on  tne  group  is  distributed  with  a  density  f (t)  : 


T  ~  fp(t) 


(2.44) 
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The  suoscript  P  indicates  the  origin  of  the  density 
(people) .  We  may  express  t  which  is  defined  as  before  as 
follows: 


t  if  no  suppressing  round 
arrives  in  (0,t)  or  t>t 

t  =  \  where  t  is  the  time  of  the 
first  suppressing  round. 
t'+  t'  if  t'<t 


(2.45) 


Remark,!1  has   the   same   distribution   as  T  .    A      similar 
derivation  as  before  leads  to  the  following  result. 


-A    t' 
s 


-X    t< 


E[x|T=t]    =   t  /   X    e      s     dt'+J      [C1  +E[x  |T=t]]X    e      S     dt< 


(2.46) 


aecall,  the  expression   A  e 
suppressing  rounds. 


-A  t 
s 


is  the  density  of  incoming 


-A   t         t  -X    t'  ^X    t 

E[T|T=tJ    =   te      S     +/      t'Xe      s     dt »   +  E[t |T=tJ (1   -   e     S) 

0  S 


,  i       -A    t  -X    t 

E[t|T=t]    =  ^-  -  ^  e      S     +  E[x|T=tJ(l   -   e      S) 
s  s 


(2.47) 


Solving    for    E[  tI  T=t],    we   find: 

E[t  |T=tJ    =  -^   (e   S     ~  1) 
Ks 

cr    replacing    A     again    by         ^s  =  A  «P(s) 


(2.48) 
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E[T|T=t]  *x7Tsy(e 


XfP(S)C 


-  X) 


(2.49) 


So  far  this  value  represents  the  expected  duration  spent  in 
the  "down  state"  given  a  particular  time  gap  t  for  a  certain 
individual.  Clearly  the  conditional  expectation  has  the 
same  fcrtn  as  the  unconditional  expectation  E[t]  (see  formula 
(2.43))  for  a  fixed  time  gap  T. 

Removing  the  condition  in  equation  (2. 49)  we  find 
the  expected  duration  time  of  suppression  based  on  the 
considered  population. 


E[t]  =  /  E[T|T=t]-f  (t)dc 
0 


,  X    P(S)t 

EtT]  ■/  r^k  (e  -D-fP(Odt 


(2.50) 


(2.51) 


E[t]    = 


XfP 


-     XfP(S)t 

<k{  e        fp(t)dt- 


XfP(S) 


(2.52) 


For  an  evaluation  of  this  expectation,  the  density  f  (t) 

P 
to  be  known. 


has 


A  further  variation  of  this   process  leads   to   the 

second   approach.    Here  the  time   gap   T  for   a   certain 

individual   varies   randomly   according   to  a    particular 

distribution.    This   approach  is  limited  to  one  individual, 
and  will  not  be  extended  to  a  group. 
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Suppose  an  individual  reacts  with  a  random  time  gap  t 


T=t, 


X"  I  t  *.  0  -J  /  •  •  • 


en   any   incoming   suppressing  round  i,  and  assume  that  this 
random  variable  T  is  distributed  with 


T  ~  fT(t) 


(2.53) 


The  subscript  I  indicates  the  origin  cf  the  density  (time). 
The  duration  x  is  defined  as  in  formula  (2.45)  before.  The 
conditional  expectation  can  be  written  as: 


-X   t         t 
s 


-X    t' 
s 


E[t|T  =  t]    =    te      5     +  /      (t»   +  E[t])    X    e      °      dt1 

os 


(2.54) 


Ihis    time     T     is    conditioned    on    the    individual's    first    chosen 
time   gap    1. 


'  -X   t  -X   t 

EfxlT^t]    =  -±-  (1  -  e      s   )   +  E[t](1-   e      s 
s 


) 


(2.55) 


Removing  the  condition  in  equation  (2.55)  we  can  express  the 
expected  duration  time  of  suppression  based  on  the 
individual's   time    gap   distribution. 


E[t]    =  /      E[t|T  =t]-f    (t)dt 
0  L 


(2.56) 


E[x]    =  / 


.  -X    t  -X    t 

~  (1  -  e     s   )   +  E[t](1  -  e     s   ) 

s 


fT(t)dt        (2.57) 
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Solving    tor    the    expected   value   E[ t  ]   and   inserting 

Xe  =  X,«P(S) 
s  £ 


(2.36) 


we   find: 


-X    P(S)t 
/      (1  -   e  )    f    (t)dt 


E[t] 


(2.58) 


XfP(S) 


[w° 

L        o 


-X   P(S)t 
(1  -   e  )    fT(t)ck 


] 


cr    expressing   in   Laplace   transform  with   s   as   an   argument: 


E[t]   = 


XfP(S) 


f(s) 


-  1 


(2.59) 


A  further  evaluation  of  this  value  requires  the  distribution 

yt)  of  i. 

comparing    toth   expectations, 

Expected      duration      time   of   suppression    based    en    the 
population 


E[t]    = 


XfP(S) 


0 


«    x  p(s)c 

e  fp(t)dt   -   1 


(2.52) 


and      expected      duration    time   of   suppression   based   on 
the    individual's   time   gap   distribution 


E[t]    = 


X-fP(S) 


f(s) 


-  1 


(2.59) 
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lie   observe   that   they  are  different.  However  if  we  suppose 

that  both  distributions  f (t)  and  f(t)   are   concentrated   at 

P  .       T 
T  =  z    (delta    function)   it   is   possible   to   reduce   both 


expressions  to  the  very  first  derivation  (2.48) 
the  time  gap  1  was  fixed. 


where   the 


The  reciprocals  of  these  expectations  approximate  the  rate 
at  which  the  foxhole  occupant  (e.  g.  member  of  a  group  of 
antitank  gunners  returns  from  the  suppressed  state  into  the 
unsuppressed  state)  returns  back  to  continue  his  mission. 

1 


Au  ~  E[t] 


where  E( t  ]  represents  any  of  the  derived  expectations. 
In  summary  the  three  rate  coefficients  developed  are: 


(2.60) 


Ak  =  Xf-P(K) 


Xs   =  WS) 


A   = 


u    E[t] 


(2.38) 
(2.37) 
(2.60) 


If  killing  as  an  additional  event  is  considered,  the  rate  of 
rise  is  computed  with  the  same  formula  (2.60)  except  that 
this  time  P„(S)  is  used  instead  of  P(S)  . 


C.   SUPPRESSION  MODEL 

The  model  which  will  be  developed  in  this  section  can  be 
set  schematically  in  the  following  framework: 
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SCHEMATIC    REPRESENTATION    OF    THE    SITUATION    TO    BE 
KCDELED 
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This  analysis  is  restricted  to  the  effects  en  the  anti-tank 
gunners...  in  this  case  the  X-forces.  It  is  behavior  of 
members  of  the  X-force  that  was  uodeled  previcisly. 

The  approach  to  formulate  the  situation  of  the  anti-tank 
gunners  is  done  by  using  the  ideas  of  Lanchester-type 
eguations   and  their  further  development  [ Tayler, 1978 , p. 20  ]. 

The  differential  equations  representing  the  model  are 
all  deterministic  in  the  sense  that  each  of  them  will  always 
yield  the  same  output  for  a  given  set  of  input  data.  Even 
though  combat  between  military  forces  is  a  complex  random 
process.  These  equations  shed  light  on  combat  dynamics  and 
may  oe  useful  in  defense  planning  studies. 

Ihe  basic  idea  is  that  artillery  forces  use  "area"-fire 
tc  suppress  or  eliminate  forces  like  anti-tank  gunners. 
"Area"  in  this  context  means  the  fact  that 

the  artillery  unit  "knows"  tne  area  in  which  to 
shoot,  but  does  not  know  the  location  of  each 
anti-tank  gunner, 

the  anti-tank  gunners  are  "invisible"  to  the 
artillery  unit. 

If  we  further  assume  homogeneous  forces  of  X,  it  is  possible 
to  set  up  differential  equations  which  model  the  rate  of 
change  of  the  X-forces: 


dX  (t) 

— rf =  -X  X   (t)  -  A.X  (t)  +  X  X  (t) 

at        s  a       k.  a       us 


(2.61) 


dX  (t) 

— f =  X  X  (t)  -  X  X  (t) 

dt_       s  a       us 


(2.62) 


Bemark:  Clearly  X  here  is  taken  from  the  4 -event-model .   The 

s 
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variables  X^t)  and  Xg(t)  are  respectively  the  number  of  the 
X-forces  which  are  either  active  in  the  foxhole  (able  to  use 
guns  and  anti-tank-weapons)  ,  jyt)  ,  or  suppressed  ,  Xs(t). 

In  this  simplified  structure,  the  two  equations  describe 
the  most  essential  factors  of  the  assumed  situation.  They 
are  mathematically  approximate  ,  because  the  solution  of  the 
two  differential  equations  will  furnish  the  number  of 
gunners  active  on  the  battlefield  and  the  number  of  gunners 
suppressed  at  any  point  in  time.  The  ^  derivation  uses  the 
Laplace  transformation: 


£{i£l)    .   „<„  -  X(0> 
^  {X(t)}   =  x(s) 


(2.63) 


ii€  perform  Laplace  transfor maticn  upon  (2.61)  and  (2.62) 
sx  (s)  -  X  (0)  =  -  (X  +  X  )  x  (s)  +  Xuxs(s) 


(2.64) 


sx    (s)    -   X    (0)    =   Xsxa(s)    -   Auxs(s) 


(2.65) 


Solving  these  two  equations  for  XJ-)  and  Xg(t)  and 
translating  them  aack  to  X&(t)  and  X^t)  gives  the  desired 
time   dependent   quantities: 


-Xa(0)    =  -   (Xs  +  Xk  +  s)    xa(s)    +  Xuxs(s) 


(2.66) 
(2.67) 


-Xs(0)    =     Xsxa(s)    -    (Xu   +  s)    xs(s) 

-Xa(0)[Xu+  s]    =   -(Xg   +  XK+   s)(X      +   s)    xa(s)    +  Xu(Xu  +   s)    xs(s) 

-Xs(0)    Xu  =   XsXuxa(s)    -    (Xu  +   s)Xu   xs(s)' 

-Xa(0)(Xu  +  s)    -   Xs(0)Xu   -   XsXuXa(s)    -    (Xs  +  Xk  +   s)(Xu    +  .)    xfl(s) 

xjo)  xu  +  xs(Q)  xu  +  sxa(0)  (2-68) 


x   (s) 
a 


(X    +  A,  +s)(X     +s)  -     XX 
s       k  u  s  u 
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multiplying  (2.  66)  and  (2.67)  with        Ag         and 

(X     +  X,    +  s)  respectively  we  receive: 

S     K 

X  (0)X   +  X  (0)(X   +  X,  )  +  sX  (0) 

xs(s)  =  (x  +x1t  +  sH\    -fs)  -x  x (2,69) 

s    k       vi         s  u 


Suppose  that  at  t=0  the  number  of  active   gunners   (able   to 

watch   and   shoot)   is   equal  to  X  and  the  number  of  gunners 

being  suppressed  (down  in  the  foxhole)  is  equal  to  0, 

for  i.e.  t=0 

X  (0)  =  x 
a      u  (2.70) 

X  (0)  =  0 
s 

the  equations  for  X  (s)  and  XJs)  can  be  rewritten: 

3.  S 

XnX  Xns 

x   (s)    =  -, °-^ + °- (2.71) 

s     +  s(X      +  A,    +  X    )    +  X,  A  s      +   s(X      +  X.    +  X    )    +  X.  X 

sku  ku  sku  ku 

Vs 
x   (s)    =  -5 ^ (2.72) 

S  s     +  s(X      +  X,    +  X    )    +  X,  X 

s  k  u  k  u 


Using   the   Laplace    correspondence: 


f    leat-ebtl  1 

<*-     )       a  -  b        I    ~  1  for     a  r  b 

v  '         s     +    (-a-b)s   +  ab 


(2.73) 


£ 


at     ,    bt   ) 
ae      -  be        ( 


a  "  b        12 

s     +   (-a-b)s  +  ab 


for     a  ^  b 
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where 


°  k  u 


and  -a  -  b  =  X     +  X.    +  X 
s  k  u 


(2.74) 


hence        a  =  ~[-(X     +  X.    +  X   )  ]  +  /(Xe  +  X.    +  X    )2  -  Ax.  X 
*  I  s  k  u  s  k  u  ku 


and 


tC-(a    +  x,   +  X  )]  _  J(x    +  x.   +  x  )2  -  AX.  X 

2sku  sku  ku 


(2.75) 


Ihe    equation    for    X     (t)     and    X    (t)     then    are: 

a  s 


X    (t)    =  \      [(X      +  a)eat   -    (X      +  b)ebt] 

a  a   -    b  u  u 


(2.76) 


S  a   -    b        s 


(2.77) 


Since  \   .  A,   and   X   have  tc  be  always  greater  or  equal  to 
s    k        u 

zero 


(2.78) 


X 
s 

>_ 

0 

Xk 

> 

0 

X 
u 

> 

0 

The   constants  a  and  b  are  always  negative  and  real  numbers. 

a  <_  0 

b  £  0 
This  leads  to  the  basic  shape  of  the  function 


(2.79) 


X  (t)  =  f(c) 
a 


X  (t)  =  f(t) 
s 


(2.80) 


shown  in  the  following  figure. 
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A  xa(t) 


x    -t 


0 


A   x  (t) 


Figure    9    -    FUNCTIONS  X   (t)    AND  X  (t) 

a  3 
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In      discussing      the      functions      X  (t)       and    X  (t)    tha    following 

a  s 

properties   can   be   seen. 


for     t  -  0 


V0)    -  a-^o"   "*u   +  *>*"°  "    <AU  +  ^'^    =  X0  (2-81) 


X    (0)    = X    [e  -   e         ]    =   0 

s  a  -   b        s 


(2.82) 


for      t   =•  °° 


X    («•)    =    lim    ~    [(x      +  a>eat-    (X,,   +  b)-ebt]    =   0        (2.83) 

a  a  -   b  u  u 


X    (»)    =    lim X    [e        -    e  =  0 

s  a  -   b        s 

t  •*■<» 


(2.84) 


for    t   =t  (i.e. max    values    for    X    (t)     and    X     (t)) 

max  a  s 


X   (t)      has   no   max   in      (0,°°),    i.e.    max   value    is   at      C   =  0 
a 


V°  =  xo 


(2.85) 


X    (t)      has   a  max  at 


£n(n/b) 
'max  b   -   a 


w  i  t  h 


(2.86) 


X_  at  bt 

X    (t        )    =  °-x    [e     "iaX  -   e     n,aX]    <   X(2.87) 

s      max  a   -    b      s  U 


Ihe  foregoing  model  assumed  a  rate  of  killing  in  addition  to 
a  suppression  rate  and  a  rising  rate.  It  is  possible  to 
simplify  this  model  by  leaving  out  the  third  rate.  It  will 
ce    shown    in    chapter    III.  C,    that    this   rate 


Ak  =  Xf-P(K) 


(2.38) 


55 


for  such  a  scenario  may  be  very  small  in  comparison  to  the 
other  two  rates. 

For  this  reason  we  can  state  that  leaving  cut  the 
killing  rate  will  not  drastically  oversimplify  the  codel, 
yet  it  will  simplify  the  ccmputat icnal  procedure  in 
obtaining  the  wanted  dependency  between  time  and  the  number 
of  gunners  active  or  suppressed  on  the  field. 

This  leads  to  the  following  set-up  of  differential 
equations : 


dX.(t) 


dt 


-A  X  (t)  +  A  X  (t) 
s  a       us 


(2.88) 


dX  (t) 


dt 


=   A  X  (t)  -  A  X  (t) 
s  a       us 


(2.89) 


where  the  total  sum  of  people  either  suppressed  or  active  on 

the  battle  field  is  equal  to  a  constant  X   (no  killing) 

o 

i.e.  X  =2  (t)+X  (t) 

o   a     s 


which   enables    us    to   rewrite    the    equations: 


dX  (t) 

— J— -  =  -A   X   (t)   +  A    (X_   -  X   (t)) 
dt  s   a  u     0  a 


(2.90) 


dX    (t) 

f-        =      X    (X.    -   X    (t))    -    A    X    (t) 
dt  s      0  s  us 


(2.91) 


Using   again    Laplace   transformation   as   earlier. 


A  X 

sX    (s)    -   X.(0)    =   -A    X    (s)    +  ~"  -    A    X    (s) 
a  0  s   a  s  a   a 


(2.92) 


A   X 

sX    (s)    -   X    (0)    =  — S_Si  -    A    x(s)    -    A    X    (s) 
s  s  s  ss  us 


(2.93) 
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x  xn  X  (0) 

Vb;     '   s[s  +   (X      +  A    )]        s  +   (A      +  X    ) 


s         u 


s  u 


x>;x0  xa(0) 

Xs(s)    =   s[s  +   (A     +  X   )]    +  s  +    (A*    +  X    ) 


s  u 


s         u 


using   the    Laplace    transform 


;£ii  «-"">!  - 


1 


j  s(s  +  a) 


.£<•-"> -ri 


s  +  a 


(2.94) 


for   both  agnations  we  find  X  (t)  and  X  (t)  fee  the  simplified 

3.  S 

model. 

The  tatt  As  has  to  be  derived  this  time  frcm  the 
probability  cf  suppression  P  (S)  without  Killing  included  as 
a    possible    event.       We    recall; 


P(S)    = 


K 


,     (2.19)       A      -  X      •    P(S)  (2.36) 

s  r 


2-iTO 


1  + 


-    E 


2a 


Hence 

X    (t)     =       A    Xn ±—    (1   -   exp[-(A      +  A    )t])    +  X    (0)    exp[-.(\    +  A,  )t] 

a  uUA+A  s  u  a  s        u 

s  u 

(2.95) 


Xg(t)     =       A_X 

for    t   =    0 


s"0   A^T-   (1  '   eXp[-(Xs  +  Vt])    +  Xs(0)    ^f-<^s  + Vt] 


X  (0)    =  xn,  X   (0)    =  0 

a  0  s 

I  he   equation   for   X  (t)    and   X  (t)    then  are: 

a  s 


(2.96) 
(2.97) 


0 


x  (O  =    -, 
a  A     +  A 

s  u 

xs(t)  =  xo  T-TT 


[A     +  A  •  exp[-(X     +  X   )t]] 
us  s  u 


(1   -    cxp[-(As  +  Au)t]) 


(2.98) 


(2.99) 
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Bemark:   For 


the  sua  of  both  assymptotic  values  of 


X(t)  and  X(t)  is  equal  to  X  . 
as  0 

The   basic  shape  of  these  two  functions  X(t)  and  X (t)  can  be 

a         s 

seen  in  the  following   Fig  10. 

If   we   compare  this  result  to  the  result  on  figure  9  we 

can  see  that  the  addition  of  killing  to   the   model   has   an 

influence   on   the   shape   of  the  functions.  So  is  e.  g.  the 

value  for  X(t)  and  X(t)  in  the  first  model   (figure   9)   for 

large   times   approximately   zero,   while   in  this  case  here 

(figure  10),  the  total  number  of  people  on  the  field  (X(t)  + 

X(t))  is  always  constant,  i.  e.  both  functions  approaches  to 
s 

a  limit  value  not  equal  to  zero. 
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A  *a(fc) 


A  x  (t) 


Figure    10    -     FUNCTIONS  X   (t)   AND  X  (t) 

a  3 
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III. 


FITTING  EXPERIMENTAL  DAIA  DESCRIBING  SUPPRESSION 


A. 


EXPERIMENTAL     DESIGN 


1 


General  Aspects 


The  mathematical  model  developed  in  the  chapter  before 
shall  now  be  supported  by  an  experiment  which  was  conducted 
by   the   US  Army  Combat  Developments  Experimentation  Command 

(CDSC)  ,  fcrt  Ord,  California. 

The  experiment  described  here  belongs  to  a  series  of 
similar  experiments,  all  dealing  with  the  objective  of 
collecting  and  analyzing  data  of  the  suppression  process. 
The  data  of  this  particular  part  of  the  CDEC  experiment  were 
based  on  the  relationship  between  suppression  and  distance. 
Ihis  analysis  tries  to  make  use  of  the  data  by  analyzing  the 
relationship  among  suppression,  distance,  and  angle.  In  a 
further  experiment  conducted  by  CDEC,  tne  objective  was  also 
to  include  the  angle  as  an  additional  variable  for  the 
suppression  effect. 

2«   Setup  And  Realisation 

The  data  were  taken  from  a  part  of  the  experiment  which  was 
executed  at  Ecrt  Hunter  Liggett,  California. 

four  foxhole-bunkers  were  constructed  which 
guaranteed  the  safety  of  the  players  as  well  as  reproduced 
the  real  scenario  as  close  as  possible.  Their  tops  were 
celow  ground  level,  and  covered  with  several  layers  of  wire 
mesh  and  steel  plates.  This  provided  overhead  protection 
from  fragmentation.  Each  bunker  was  equipped  with  a  mirror 
system  which  allowed  the  player  to  view  events  in   front   of 
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his  position. 

A  pop-up  silhouette  was  installed  forward  of  each 
player  position.  The  player  was  able  tc  control  the 
silhouette  as  well  as  the  cover  of  the  mirror  system,  i.e. 
when  the  cover  was  opened  (allowing  the  player  to  view 
down-range),  the  silhouette  was  in  exposed  posture.  He  was 
asked  to  respond  as  if  he  would  be  in  the  position  of  the 
pop-up  silhouette. 

Each  bunker  was  connected  to  controll  bunkers  by 
communication  and  instrumentation  wires  and  power  bunkers 
for  data  recording  and  supply.  In  the  forefield  of  the 
bunkers,  different  types  of  projectiles  or  equivalent 
charges  (31  mm, 105  mm, and  155  mm  among  others)  were  placed 
and  randomly  detonated  one  at  a  time  with  the  time  between 
detonations  randomly  selected  from  three  possibilities  of 
ten,  fifteen,  and  twenty  seconds.  The  figure  on  the  next 
page  shows  the  schematic  setup  of  the  rounds,  the  location 
of  the  bunkers  and  the  angle  of  sight  for  each  foxhole.  It 
can  be  seen  that  the  explosions  were  visible  to  some  players 
but  not  to  others.  Since  each  type  of  ammunition  has  a 
different  lethal  radius,  it  was  necessary  to  have  different 
range  configurations  for  each  type.  The  aspect  angle  and 
the  miss  distances  from  the  foxholes  to  the  different  points 
of  detonation  are  summarized  in  appendix  A  for  each 
ammunition  separately.  Since  all  of  the  projectiles  used  in 
this  part  of  the  experiment  were  statically  detonated,  it 
a as  not  possible  to  model  the  kinetic  contribution  to  the 
terminal  effects.  In  order  to  keep  the  fragmentation  pattern 
as  close  as  possible  to  those  of  incoming  rounds,  the 
projectiles  were  supported  with  sandbags  [  5 uppex, p. A- 1 2  ]. 
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Figure  11  -  SCHEMATIC  RANGE  FOR  81/105/155  mm 
SUPPRESSION  EXPERIMENT  (SUPPEX  II) 
CDEC 
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The  players  were  divided  into  two  four-man  teams* 
For  each  trial,  members  of  one  team  occupied  individual 
foxholes  and  provided  the  performance  data  for  that  trial. 
The  mission  of  the  players  was  to  track  moving  target  tanks 
by  operating  the  periscopes.  This  mission  was  interrupted  by 
the  players  responses  to  detonations  in  case  he  was 
suppressed  (change  of  his  state  to  "down") .  The  state 
change  and  the  pericd  of  suppression  were  automatically 
recorded . 

3.   fiesentaticn  Of  The  Data 

Appendix  A  summarizes  the  data  which  are  the  basis  for  the 
succeeding  data  analysis  and  serve  tc  evaluate  the 
parameters  and  rates  for  the  model. 


E.    DATA  ANALYSIS 


1  •   Parameters  Of  The  Model 

The  data  presented  in  the  section  above  are  separated   into 
three  different  sets: 

Set  one  consists  of  data  related  to  the  conditional 
freguency  of  suppression  P(5/9,r)  and  the  time  of 
suppression  when  rounds  of  caliber  81  mm  were  fired. 
Specifically,  if  n  trials  were  made  under  conditions 
(r,9),  the  number  s.cf  suppression  was  tabulated. 
Then  s./n  is  an  estimate  of  P  (S/r.  8.,). 

Sets  two  and  three  provide  the  same  data  except  that 
they  are  related  to  105  mmm  rounds  and  155  mm 
rounds. 

The  main  part  of  this  analysis  is  to  make  use  of  the  data  in 
order  to  estimate  the  parameters  of  the  model   described   in 
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chapter  II,  and  to  see  how  well  the  real  world  situation  can 
te  described  by  the  model  developed  earlier. 

The  model  is  fitted  (numerical  values  of  the 
parameters  are  determined)  by  the  method  of  least  squares. 
[De  Groot,  p. 527-538].  With  the  basic  model  being  of  the 
form: 


Yi  =  6lUi  +  62Vi  +  Erri 


(3.1) 


And   the   data   are  analyzed  on  the  following  two  different 
assumptions;  the  second  is  certainly  the  most  realistic: 


Homoscedas ticity   of   the   data,   i.   e.   the 
variance  of  Err  is  assumed  to  ba  constant. 


error 


Heteroscedasticity  of  the  data,  i.  e.  we  will  assume 
that  each  error  term  Err  is  distributed  with 
variance  a2  , where  the  variance  is  not  constant  over 
observations.  Errors  are  also  assumed  tc  be 
independent.  If  convenient,  error  terms  will  be 
assuned  to  be  approximately  normally  distributed. 

These  assumptions  imply  different  regression 
methods.  For  the  first  assumption,  the  result  of  unweighted 
single  step  and  iterative  regression  methods  will  be 
presented.  On  the  assumption  of  heteroscedasticity,  two 
special  iterative  methods  will  be  used.  When  we  allcw  for 
heteroscedasticity,  ordinary  least-squares  estimation  places 
the  same  weight  on  the  observations  which  have  small  error 
variances  as  on  those  with  large  error  variances.  By 
applying  a  weighting  regression,  it  is  possible  to  adjust 
for  the  heteroscedasticity.  So  the  two  announced  methods  for 
the  second  assumption  are  iteratively  weighted  least  sguares 
regression  methods. 
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In  order  to  prepare  the  formula  of  the  model  fcr  the 

regression,  we  apply  a  log-transformation   to   the   equation 

(2.1)    in   chapter   II   B.1.   This  is  a  convenience,  for  it 

transforms  the  problem  to  one  of  linear  fitting.   We   obtain 

the  following  result: 


P(S/9,r)  =  exp[-  -  r  (1  -  e  cos  9)  1 


(2.1) 


Jin  P(S/8,r) 


1   2  ^  e   2 
-  -  r  +  -  r  cos  0 


(3.2) 


Using  the  notation 


y  =  In  P  (S/9,r) 


u  =  r 


(3.3) 


v  =r-cos9 

the   equation   can    be   applied    to      each      datapoint      i      and      the 
aquation    can    be    rewritten -as: 


yi  =  "  KMi  +  KVi 


(3.4) 


Ihe   model   car   then   oe    expressed   as   follows; 


Yi  =  8lUi  +  B2vi  +  Erri 


(3.1) 


where  the  unknown  parameters  are:   8-  =  -  —      ^2  =  K   (3-5) 


65 


and  the  random  variable 

Y  =  £n  P(S/6i,rjL)  (3-6> 

Ihe  "hat"  on  the  letter  P  represents  an  estimate  of  the 
probability  of  suppression  at  (r,9).  The  fitted  regression 
plane  for  this  function  has  the  form : 


{jU  +  B2v  (3.7) 


The  log-transf crmation  of  the  conditional  frequency 
cf  suppression  P(S/9,r)  causes  difficulties  because  of 
experimental  results,  which  lead  to  an  observed  frequency  of 
zero  suppression.  This  fact  influences  the  estimate 
E(S/8,r),  that  is  used. 

Having  specified  the  model  it  is  possible  to  apply 
the  following  different  regression  methods  en  the  assumption 
of  homoscedasticity: 

In   other   words,   suppose   ni  observations   were   made   at 

experimental  condxtions  (r.,  9.)  and  of  these  s  were  successes. 

11  i 

Then    first   consider   the  estimate 

Y     =   Jin  P(S/r   ,8    )    =   Zn(s./n.)  (3-8) 

where   Sj^  is  the  number  of  suppressions,  and  n^  is  the  number 
of  exposures  to  suppression,  under  condition  i. 

Consequently,   the   transformed  response  to  conditions  (r-,  84.) 
may  be  In  (9/ni)  =  -  °°    formally,   causing   embarassment.   We 
will   subsequently   suggest  alternative  ways  of  dealing  with 
this  problem. 
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Hethod  1:   (zero-probabilities  are  omitted) 


With  this  method  ,  the  datapoints  which  had  a 
frequency  of  suppression  P(S/r,6)=s7n  =0  are 
deleted.  Consequently  the  followinq  data  points  of 
the  oriqinal  sample  of  appendix  A  as  shown  in'  figure 
Fig  12  were  net  considered. 
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Caliber 


No. 


P(S/G,r) 


9 


81  mm 


29 
30 
31 
32 
34 
35 
36 


0 
0 
0 
0 
0 
0 
0 


41 

36 

60.5 

46 

60.5 

46 

36 


117 
127 

95 
108 

95 
108 
127 


105  mm 


15 
19 
23 
24 
28 
29 
31 


0 
0 
0 
0 
0 
0 
0 


17 
36 
46 
46 
56 
56 
71 


104 
158 
135 
135 
108 
108 
97 


155  mm 


3 
15 
17 
18 


0 
0 
0 
0 


0 
14 
23 
14 


200 
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Figure  12  -   DELETED  DATA  POINTS 
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Method  2  and  3  (Clustering  of  data  points) 


One   or  more  zero  data  points  i  ,  i.e.   P(S/e,r)  =  0 

i=1,2,...m   will  be  clustered  with   one   other   data 

pcint  j  which  has  a  nonzero  conditional  frequency  of 

suppression,  i.e.   P(S/9,r)  *    0  and  which   possesses 

the   closest  distance  r  and  angle  9  tc  the  zerc  data 

J  J 

pcint  (s)  . 

Ey  clustering  (m+1)  data  points  together  and  taking 
the  mean  of  the  frequencies  ,  9  and  r  within  each 
cluster,  we  are  able  tc  replace  the  data  points  of 
the  cluster  sets.  This  procedure  is  inevitably 
somewhat  arbitrary. 


P(s/e,r)  =  ^-  I    P^S/e.r) 


1=1 


1     m+1  j  m+1  (3-9) 

=  m+I  I     9i         r  =  m+1  Jx  ri 


The  bars  above  the  symbols  represent  the  cluster 
average.  With  these  averages  two  different 
replacement  procedures  can  be  applied: 

Each  cluster  will  be  replaced  by  its  cluster 
average.  This  reduces  the  data  sample. 

Each  cluster  point  will  be  replaced  by  its  cluster 
average.  Consequently  the  numoer  of  data  pcints 
stays  unchanged. 

figure  13  lists  the  replaced  data  points  of  the   samples  of 
appendix  A. 


69 


O 

X. 

4-1 

oj 


(0 

4) 
3 
t-l 
CO 

> 
55 


u 

CD 
CO 

P_ 


r-~  r»  f»» 

CM    CM    CM 


1/1    1/1    m    in 
CM   CM    CM    CM 


o  o  o  o 


m  m  u*i 

CTi  ON  <r> 


o  o 


co  en 


o  o  o 
co  co  co 


cm  cm  cm  cm 

oooo 


o  o 

o  o 

CM    CM 


r~-  f^  r-~      10  m 

O    O    O        CM    CM 
CM    CM    CN        i — I    i — I 


VO    vO   v£> 
CO   CO   CO 


OOOO 
-3-    <•    -J   <T 


O    O    O 
\0  \0   \D 


KO    M3 
CO    CO 


r~  t^  t~- 

<r  <r  <r 


<r  -J-  <r  <f 

vO  vO   vO  \o 


O    O        CM    CM    CM 


CO    CO 
CM    CM 


CO    CO   CO        CO    CO    CO    CO 

•^  <r  <r      co  co  co  co 
o  o  o     o  o  o  o 


r-~  r-»  r-. 

in  uo 

m  1/1 

\0   vO  \o 

vD  vO 

\0  vO  MD 

CO  CO  CO  CO 

01  on 

■ — (  • — 1  ■ — 4 

O  O 

i-l  i-l 

o  o  o 

o  o  o  o 

o  o 

co  co  co 


o  o 


ooo  oooo  ooo  oo  oo  ooo  oooo  oo  ooo  oo 


1/1 

CM 


UO 

ON 


o 


co 


o 

CO 


CM 

o 


o 
o 

CM 


o 

CM 


CM 


o 

u 
X 


CO 


O 


o 


CO 


-3" 


CD 

co 
On 


co 
o 


co 
co 
O 


m 

M3 

O 


o 


CO 

o 


On 
O 


CO 
CM 


O 


r-»  r-.  r-» 

CM    CM    CM 


co  co  r--  oo 

O  O   i-l  o 


1/1  1/1  1/1 

ci  ci  oi 


o  o 


CO   O 

UO     r— I 


i/1  i/1  o 
CO    CO    CM 


CO  CO  r-«  r- 

o  o  oi  oi 


o  o 
o  o 

.si    CM 


CM  OI  On 
OOO 
CM    CM    CM 


in    ul 
CM    CN 


CO 
V 

3 


O   >0  vO 
CO  CO  CO 


vo  ^o  i— i  r^ 

-T    <■    -3-    CM 


in  in  i/i 


OOO      r--  r-~ 

\D    \0    O        i-H    i— * 


*n  r-~ 

\D  \D   ON 

vO  \0     r—l  f— 1 

co  co 

-3"  <"  <T 

in  m  r-»  r-~ 

O  O      r-~  -3-  <r 


CO     CO 
CM    CM 


U 

n 

CD 

CO 


CO 


OI       00 

r-H         CO 


rH  UO        i-l  CM  rH  >-( 

OOO   OOOO   OOO   OO   OO   OOO   OCT.OO   OO   OOO   OO 


o 


co  o  no  in  m  en  •*  <j  h  co  <r  m  oo  co  <r  m  oo  oi  o  •— i  coo  oimco  <r  r— 

COCOCO        COCOCMCM        COCOCM        .-HrH        H   M        CM    CM    CN        CSICMCOCO  iHr-l        r-H    r-H 


U 
X> 

•w 

(0 

c_> 


m 

o 


m 
m 


Figure  13  -   CHANGED  DATA  POINTS 
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Method  4  (data  transformation  by  Cox) 


This  method  is  analogous  to  one  suggested  by   D.   R. 
Ccx  [Cox,  p. 33].   The  derivation  for  the  variable  yi 
shewn  on  these  pages  can  be  applied   correspondingly 
in  the  following  manner: 


Y.  =  In 
~i 


S  +  a 


(3.10) 


where   the   random   number   Si  is   the    number   of 
suppressions  (successes,  given  ni trials) . 


Ihe  constant  a,  which  represents  an  unfciasing 
adjustment,  is  derived  in  appendix  E.  It  is  taken  to 
te: 


a  =  i  (1  -  P(S/8,r)) 


(3.11) 


This   is   suggested  by  an  auxiliary  analysis  similar 
to  that  of  Cox. 

Hence  the  transformation  results  in  the  formula: 


Xi 


=  £n 


S±  +  Y  (1  -  Pi(S/0,r)) 


(3.12) 


This  formula  enables  one  to  include  data  points  that 
involve  zero  observed  frequencies  cf  suppression; 
the  embarassment  of  taking  the  logarithm  of  zero  is 
no  longer  present. 
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Method  5  (unweighted  iterative  regression) 


As   an   extention  of  the  method  discussed  last,  this 
method  replaces  the  value  of  P(S/6,r)  in  the  formula 


S,  •+  \   (1  -  P.  (s/e,r)) 

Y\    =    (  — - }  (3.13) 

~i    \        n. 


After   each   iteration   by   an   estimate   of  the 

probability  of  suppression  P  (S/9,r)  ,  the   number  of 

iterations   was   determined   by   the   appearance  of 
convergence  of  the  parameters  K  and  e. 


The  two  following  methods  are  weigatad  regressions  which  are 
necessitated  ty  the  earlier  described  ae tercscedast icity. 


Method  6  (weighted  iterative  regression  Var(Zrr)~r  ) 


This  method  is  suggested  by  the  fact  that  the 
observed  variability  in  residuals  increases  with  r. 
[Pindyck,  p.  1C0]  Since  the  error  variance  is  not 
known,  it  is  reasonable  to  assume  the  existence  of  a 
simple  relationship  between  the  error  variances 
Var  (Zrri)  and  the  values  of  one  of  tne  explanatory 
variables  in  the  regression  model.  In  this 
analysis,  the  distance  from  tne  foxhole  to  the 
explosion  was  chosen  as  the  important  explanatory 
variable.   By  using 
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Var[Err   ]  ~  r 


(3.14) 


the    resulting    regression   formula   is: 


(Var[Erril) 


1/2 


(Var[Err±]) 
Err. 


1/2  +  B2 


(VarfErrJ) 


1/2 


(Var[Erri]) 


1/2 


(3.15) 


This      equation      can      be      reduced      to      the      following 
regression   function: 


y'.   =  3,   +  B-v! 
i  1  2    i 


(3.16) 


where 


YI  = 


in 


S.   +  j(l  -  P.(S/d,r)) 


(3.17) 


•i-5 


Afterwards,  the  iteration  procedure  is  equivalent  to 
method  5.  This  method  tends  to  make  the  variance  of 
the  residuals  around  the  fitted  line  of  more  nearly 
constant  variance;  estimates  of  the  model  parameters 
should  be  thereby  improved. 


Method  7  (weighted  iterative  regression) 
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Using  the  error   variance  derived   in   appendix   C, 
which  has  the  form; 


VarUn  f,] 


1  -  Pi(S/9,r) 


(3.18) 


where 

f,  =  —  (3.19) 

1   ni 

the  regression  formula  will  be  modified  as   follows: 


=  3, 


(VarUn   f±])1/2  1    (VarUn   f±])1/2  2    (Var[£n   f   )]1/2 


(3.20) 


where  y  is  again: 


y.  =  In 


S.  +  y  (1  -  P.(S/9,r)) 


(3.12) 


as  derived  in  method  U .    The  resulting   equation   has 
the  form:  /         u± 

Sl  7,       r777,     s\l/2 


In 


S.  +  ±   (1  -  P.(S/8,r)) 

i   2 i 

n. 


1  -  Pi(S/9,r) 
iK-  P^S/e.r). 


1/2 


1  -  PjL(S/9,r) 
n.  •  P.(S/8,r). 


+  3. 


'l  -  P.(S/9,r) 

>n.  •  P.(S/8,r) 
x   i 


(3.21) 
T72-+  **! 


The  iterative  procedure  in  this  method  is  performed 
by  using  the  estimate  of  the  conditional  probability 
of  suppression  as  input  for  each  succeeding 
iteration.  The  number  of  iterations  was  determined 
as  in  method  5  and  6  by  the  apparent  convergence   of 
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the  estimators  of  the  parameters  K  and  e  . 


In  order  tc  perform  the  multiple  linear  regression,  we 
switch  over  to  matrix  notation,  where  we  can  write  the 
normal  equations  in  the  following  form: 


T       T 
Z  -Z-B  =  Z  «Y 


(3.22) 


The  matrix  Z  is  the  design  matrix  which  consists  of  vector 
0,  the  square  of  the  distances  and  V,  the  product  of  the 
cos8  and  the  square  of  the  distances.  The  variables  y  and  S 
represent  a  (sample  size  lxl) and  a  (1 x2)  vector  respectively, 
for  the  experiments  analyzed,  the  entries  of  the  design 
matrix  were  obtained  by  solving  the  normal  equations  for  8. 
He  find: 


T    -IT 
=  (Z-Z)  -Z-Y 


(3.23) 


Remark:  Capital  letters  used  for  matrix  notation. 

Ihe   parameters   S   and  e   can   be   evaluated   by  using  tae 
equations : 


c  = 


B2-K 


(3.24) 


As   a  modification  of  the  suppression  function  (2.1) 
in  chapter  II.  B.  1  the  function  (2.4)  in  the   same   chapter 
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was   considered. 


H 


P(S/0,r)    -  L«exp     --r    (1  -   e   cos  0) 


>] 


(2.4) 


Ihe  inclusion  of  the  parameter  L  allows  for  the  possibility 
that  suppression  may  not  occur  for  shots  that  fall  in  very 
close  proximity  to  the  foxhole. 

Method  8  (siiiilar  to  method  4  but  using  model  above) 


In  this  method  the   regression   procedure   and   the 
method  4  were  applied  to: 


Y,  =  S.u.  +  60v.  +  S  +  Err, 
i    1  l    2  i    3      i 


(3.25) 


where 


y .    =   In. 
l 


S1  +  j   (1  -  Pi(S/0,r)) 


u.  =  r. 
l    i 


v  =  r.  cos  0 


(3.26) 


Sl '  "  "  K 


2   K 


63  =  e 


The  ccrmal  equations  have  tne  same  form  as  before: 

(3.23) 


T    -1   T 
8  =  (Z-Z)   -Z-Y 
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but  1,  the  design  matrix  consists  this  time  of  the 
vectors   U,V,  and  a  column  vector  consisting  only  of 

A 

ones.  The  resulting  8  vector  consists  now  of  3 
elements,  from  which  it  is  possible  to  evaluate  the 
different  parameters  K,  s  ,  and  L. 

The  calculation  of  the  parameters  for  all  methods  was 
performed  fcr  each  type  of  ammunition  separately  by  applying 
the  computer  language  APL.  The  programs  which  where  used 
for  this  regression  are  listed  in  appendix  I.  The  appendixes 
D  through  G  show  the  analysis  of  the  residuals  after  the 
regression  cf  the  original  data  of  appendix  A. 

Ihe  analysis  for  each  ammunition  type  was  performed 
according  to  the  same  pattern,  and  includes  the  following 
steps. 

For  each  earlier  developed  method  (1  tarough  7)  : 

Determination  of  BETA(1)  and  BETA  (2) 

Determination  of  K 

Determination  of  e 

Plot  of  the  residuals  as  a  function  cf  y  as  defined 
in  each  method  earlier.  The  APL  function  which  was 
used  has  the  name  SCAT  and  belongs  to  the  library 
package  OA  2  3660  (available  at  8.  R.  Church 
Computer  Center)  . 

Plot  of  the  residuals  of  the  regression  as  a  single 
array  with  the  function  BOXPLOT  of  the  same  library 
package.  The  plot  characterizes  the  guartiles,  the 
interquartile  distance,  the  median,  data  points 
inside  and  outside  the  1  and  1.5  interquartile 
distance  and  outliers  [McNeil,  p. 13  and  71,72]. 

Numerical  values  of  the  residuals. 
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Histograms  for  the  residuals  of  the  regression  which 
show  the  relative  frequency  and  statistical 
features. 

For  method  8,  which  is  an  analysis  based  on  different 
assumptions  for  the  probability  cf  suppression  P(S/9,r),  the 
same  analysis  was  performed,  with  the  exception  of  the 
histograms  fcr  residuals. 

The  methods  1-7  describe  the  attempt,  to  master  the 
problem  of  zero-probabilities  in  the  original  data  set  as 
well  as  tha  problem  of  a  possible  heteroscedastisity . 

For  the  later  methods,  in  which  irore  appropriate 
statistical  tools  were  used,  the  systematic  structure  of  tne 
residuals  seems  to  disappear  up  to  a  certain  point.  The 
residuals  concentrate  themselves  more  and  more  symmetrically 
around  their  mean  and  median  (compare  particular  the 
boxplots  and  the  histograms  in  the  appendixes  D,  E,  F,  and 
G)  . 

The  relative  large  remaining  range  of  the  residuals 
is  determined  by  single  outliers.  The  analysis  showed 
differences  for  different  kinds  of  ammunition.  Among  the 
three  ammunitions  considered,  the  analysis  cf  81  mm  showed 
tne  smallest  spread  for  the  residuals,  for  almost  all 
regression  methods. 

The  appendix  H  shows  for  the  iterative  regression 
methods  5,6,  and  7  the  development  of  K  and  e  for  the 
different  kinds  of  ammunition.  In  each  method  ,  a 
convergence  of  the  values  K  and  e  with  increasing 
iteration  can  be  observed.  The  starting  value  for  K  is  in 
method  5  smaller  and  in  method  6  and  7  larger  than  its 
corresponding  value  after  convergence  is  obtained.  This  is 
true  for  all  types  of  ammunition.  An  equivalent  observation 
can  be  made  for  the     e -values. 

A  graph  for  the  iterating  e  for  method  5,  155  mm 
could  not  te  plotted  because  of  the  very  small  change  cf  the 
£   -values  along  each  iteration. 
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The  third  page  of  appendix  [1  displays  the  numerical 
values  for  K  and  e  produced  in  each  method  for  the  three 
different  types  of  ammunition. 

The  analysis  of  the  parameters  K  and  £  showed  the 
following  result; 

for  the  iterative  methods  (5,6,7)  the  K  and 
e-rvalues  approached  with  increasing  iteration  a  limit  value 
(see  page  one  and  two  of  appendix  H) ;  this  occured  after 
about  the  sixth  iteration.  The  final  values  for  these 
aethcds  are  also  presented  in  the  tableau  of  appendix  H 
(third  page) .  The  plots  for  K  and  e  show  different  shapes 
for  different  methods  and  partly  also  for  different  types  of 
ammunitions  .  The  approach  to  the  final  value  may  occur  from 
a  relative  snail  or  a  relative  high  value.  Among  the  three 
ammunition  types  considered,  the  analysis  of  81  mm  showed 
tne  smallest  spread  for  the  residuals,  for  almost  all 
regression  methods.  The  scale  parameter  K  influences  the 
probability  of  suppression  P(S/9,r)  as  stated  under  formula 
(2.1)   in  chapter  II.B.1.  in  the  following  manner: 

An  increasing  value  of  K  leads  to  an  increase  of  the 
probablility  of  suppression  E(S/9,r)  .It  is  reasonable  to 
assume  that  tne  suppression  effect  increases  with  the 
increase  of  the  calicer.  This  behavior  was  confirmed  by  the 
data  analysis;  except  for  method  one  and  three  where  an 
inversion  could  be  ooserved  between  the  K-values  for  S 1  mm, 
105  mm,  and  155  mm  and  for  81  mm  and  105  mm  respectively. 

This  distortion  results  from  the  fact  that  in  irethod 
one  all  data  points  with  probabilities  of  suppression  equal 
to  zero  are  disregarded  and  in  method  three  the  cluster 
procedure  emphasizes  the  average  values  produced  by  the 
clustering.   He  prefer  the  latter  methods. 

Contrary  to  these  observations  on  the  K-  values,  a 
general  trend  for  the  £ -values  can  not  be  related  to 
different  methods  or  different  kinds  of  ammunition.  The 
most   reliable  value  of  K  and  e  seems  to  be  found  by  method 
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7.  This  method  considers  the  different  variabilities  for 
different  sets  of  distance  and  angle. 

In  the  course  of  this  paper,  the  value  of  K  and  £ 
derived  in  method  7  are  taken  as  input  for  developing  the 
rates  for  the  model  of  section  II. B. 

Ihe  third  parameter  L  developed  by  method  8  is  valid 
only  up  to  a  probability  of  suppression  evaluation  equal  to 
1.0.  Since  for  31  mm  and  155  mm  L  is  larger  than  one,  the 
function  for  the  model  should  be  decomposed.  This 
alternative  approach  to  model  suppression  is  beyond  the 
scope  of  this  paper  and  will  not  be  considered  further  at 
this  point,  although  it  is  certainly  a  topic  worthy  of 
further  research. 

As  a  possible  further  step,  confidence  intervals  for 
K  and  e  could  be  developed. 

2-   Suppression  Time 

The  suppression  time  data,  which  were  collected  in 
seconds,  represent  the  time  for  which  an  individual  remains 
suppressed  as  a  reaction  of  a  singla  round.  During  this 
time,  the  individual  was  unable  to  carry  out  his  mission;  he 
is  in  hiding  in  an  effort  to  survive.  According  to  the 
setup  of  the  experiment,  it  was  possible  for  the  suppressee 
to  react  to  detonations,  which  he  was  able  to  observe  and  to 
hear  or  to  hear  only  (visual/auditory  and  auditory 
perceptors)  .  This  is  consistent  with  tha  model  design  in 
chapter  II. 

This  analysis  is  an  attempt  to  compare  the 
suppressicn  interval  gained  from  the  experiment  Kith  a 
certain  distribution  whose  parameters  are  estimates  of  the 
data.  The  GAMMA  distribution  was  selected,  because  the  range 
cf  the  random  variable  X,  representing  the  time,  is  limited 
below  by  zero  and  by  +  °°   above. 

The  data  of  suppression  interval  for  the  calibers  81 
mm,  105  mm,  and  155  mm  are  plotted  in  histograms  in  appendix 
J. 
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The  estimates  A  and  r  for  the  Gamma  distribution 
are  derived  from  the  expectation  and  the  variance  of  the 
data. 


X  - 


gjx] 

Var[x] 


i   n 

n  j    i   i 


r  =  X  •    E(x]  = 


-.2 


i=l 


n     2 

(  I     V 

i=l 

n      _  2 
n  J,   (x.-x) 

i=l   1 


(3.27) 


fcith   these    parameters,    the    Gamma    density   can    ce    computed      in 
the   considered    range,    by   the    formula: 


r(r) 


r-1     -Xx 
x  e 


0   <   x   <  « 


f(x/X,r)   = 

0  otherwise 

cor   comparison,    the    related    density    was   computed    by 

„  .    .         frequency    in   interval 

f  (x)  =  - , ; 

N  •  interval 


(3.28) 


(3.29) 


within  the  same  range,  where  the  frequency  in  a  particular 
interval  and  the  interval  itself  is  taken  from  the 
Histograms  of  appendix  J,  and  the  constant  N  is  the  sample 
size  of  the  measured  suppression  intervals. 

The  numerical  values  for  the  estimates  X  and  r  and 
the  values  for  the  sample  size  M,  the  interval  Ax,  and  the 
range  are  displayed  in  the  following  table. 


Caliber 

X 

r 

N 

INTERVAL 

Ax 

RANGE 

81  mm 

0.686 

2.174 

445 

0.5 

0  to  16.15 

105  mm 

0.468 

1.974 

348 

0.4 

0  to  13.7 

155  mm 

0.587 

2.591 

101 

0.5 

0  to  16.75 

81 


Ihe  following  Fig  14  compares  the  Gamma  density  (3.28) 
with  the  height  of  the  related  estimated  probability  cf  the 
data  (3.29).  The  sketches  show  that  the  Gamma  density  for 
each  ammunition  type  respectively  underestimates  the  related 
freguency  of  the  data  because  of  its  long  tails.  In  fact, 
any  distribution  which  meets  the  requirements  above  could  be 
taken  for  a  comparison,  although  the  Gamma  distribution 
seems  to  te  a  particularly  good  choice  for  the  81  mm  caliber 
ammunition. 
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RANGE   OF   X:    0       20 
RAN  OF.    OF   Y:    0       0.3  5 


frequency   in   interval 
H* interval 


--OO--  —  o  —  o  o  —  o—  oo  —  —  —  o  o  —  o  —  o  —  o  o  —  o  —  — 


105  mm:   A  =  0.468,   r  =  1.974,   N  =  348,   interval  =  0.4 


(r-1)    -Xx 
x      •  e 


frequency  in  interval 
M« interval 


freguoncy  in  interval 
W« interval 


• |-4 

-oo-o—  oo—  o—  o-o  ±  —  ± 


Figure  14  -   COMPARISON  GAMMA  DISTRIBUTION  AND  DATA 
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In  order  tc  intensify  the  comparison,  as  a  next  step 
the  skewness  and  the  kurtosis  of  the  data  and  of  the  Gamma 
distribution  with  the  estimates  above  will  be  contrasted. 

These  characteristics  express  the  symmetry  of  the 
distribution  about  its  point  of  central  tendency  and  the 
relative  concentration  of  cases  at  the  center  and  along  the 
tails  of  the  distribution. 

For  the  derivation  of  the  skewness  and  the  kurtosis 
for  a  Gamma  distribution,  we  need  the  first  four  moments 
about  the  origin: 


vi  =  E[x]  =  X 


r  „r  2,    f   2  -Xx  (Xx) 


r-1 


p   =  E[x  ]  =  /   x  e 
Z  0 


T(r) 


By   changing    the   variable  of    integration 
z  =  Xx 
dz  =  X  dx 

CO        2 
r,r     2,  !        Z  _Z        r~l  J  ^  L     A~ 

p2   -  E[x   ]    =  J     ^e        z  ^X-^dz 


we    find 


X     dx 
(3.30) 

(3.31) 

(3.32) 


r  (r  +  21  r      t+2-i    -z l dz 

z  e        r(r  +  2)    dZ 


J 


x  r(r)     0 

but    the    above    integral    is    equal    to   one,    hence 
U2  -  i   (r  +  1). 


(3.33) 


In  a  similar  fashion,  the  third  and  the  fourth  moment  can  de 
derived. 


y,  =  ■—  (r  +  2)(r  +  l)r 
X 

U,    =  ^-  (r  +  3)(r  +  2)(r  +  l)r 

4   x4 

Converting   these   moments  to  moments  about  the  mean 
using  binomial  expansion: 

U3  =  u3  -  3u2  E[x]  +  2(E(x])3 

u,  =  y,  -   4y.  E[x]  +  6y„(E[x])2-3(E(x]) 
4    A     3  * 


(3.34) 

by 
(3.35) 
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we  can  compute  the  skewness  and  the  kurtosis 


M/, 


'3    3 


l4    A 


(3.36) 


For  a  Gamma  distribution  we  get  the  following  formulas'. 


r.   .„    (r+2)(r+l)r     (r+l)r   r      J 

3_l     x3   "     x2    *+/V 


(rA2)3/2 


°3  = 


/r" 


and  by  similar  operation: 
3r  +  6 

Using  the  above  estimates  for  X  and  r  we  get: 


(3.37) 


a„  = 


^ 


a4  = 


3r  +  6 


(3.38) 


The   formulas  for  the  skewness  and  kurtosis  derived  from  the 

data  itself  are: 

N 


1  v ■     (  -v3 

m  2.  0 


a3  = 


(x   -  x) 

i=l 


N 


-x2N3/2 


(£  I   (x   -  x)2) 
i=l 


1  N 


-s4 


(3.39) 


a4  = 


i=l 


N 

I 
i=l 


-v2N2 


(f  I      (x  -  x)2) 


The  resulting  numerical  values  for  the  skewness  and  kurtosis 
derived  fcr  the  Gamma  distribution  and  from  the  data  are 
displayed  in  the  following  table: 
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81  mm 

105  mm 

155  mm 

Skewness 

GAMMA 

Q3 

1.356 

1.423 

1.242 

Skewness 
DATA 

S3 

2.156 

1.032 

2.103 

Kurtosis 

GAMMA 

a4 

5.759 

6.039 

5.315 

Kurtosis 
DATA 

S4 

8.004 

1.309 

6.750 

In  discussing  the  numerical  features  we  can  see  that  the 
skewness  of  the  Gamma  distribution  and  of  the  data  itself 
shows  an  asymmetric  right-skewness  (positive  values). 
Although  the  values  differ  considerably  .  So  is  e.  g.  the 
skewness  of  the  data  for  155  mm  almost  twice  as  big  as  the 
skewness  based  on  the  Gamma  distribution. 

The  same  can  be  said  for  the  kurtosis.  There  is 
significant  difference  between  their  numerical  values. 

This  analysis  supports  the  fact  that  the  Gamma 
distribution  can  only  be  a  rough  fit  to  the  data  given,  and 
as  already  earlier  expressed,  a  fit  of  another  distribution 
probably  would  have  been  as  successful  as  this  one. 


C. 


VEBIPICflllON  OF  THE  MODEL 


In  order  to  supply  general  features  as  input  factors  for 
the  decision  process  for  military  leadership,  the  purpose  of 
this  paragraph  is  to  compute  numerical  values  for  the 
probabilities   and  rates  developed  in  chapter  II.B.1.  and  2. 

For  the  computation  of  numerical  values,  it  is  necessary 
to  make  the  following  reasonable  assumptions: 

It  is  assumed  that  an  artillery  unit  consist  of  six 
weapons  all  either  105  mm  or  155  mm,  or  that  three 
81  mi  launchers  are  combined  to  a  mission  unit. 
Hence  the  fire  rate  ^f  for  such  units  will  be 
concluded  to  be  [  Wiener, p. 189, 21 1 ,2  13  ]: 

81  nir  unit  with  50,  55,  and  60  rcunds/min 

105  in  unit  with  28,  32,  and  36  rounds/min 


86 


155  mm  unit  with  18,21,  and  24  rounds/min. 

The  standard  deviation  a  for  the  density  for 
artillery  hits  will  be  considered  to  be  [  FH  6-161-1, 
p. 53]: 

25,50,and  35  m  for  81  on  unit 

r 
20,25,and  30  m  for  105  ram  unit 

30, 35, and  40  ra  for  155  mm  unit. 

The  numbers  are  taken  from  the  field  manual  6-161-1 
and  out  of  working  papers  of  CDEC.  Ihey  are  rounded 
fcr  convenience. 

The  parameter  H  of  the  probability  of  kill  is  chcsen 
completely  arbitrarily  with  100  for  81  mm,  150  for 
105  mm,  and  200  for  155  mm  and  has  nothing  to  do 
with  experimentally  observed  values  fcr  the  weapons 
here  in  question. 

The  parameter  estimates  K  and  e  are  taken  from 
appendix  H  (tableau)  with  the  values: 


81  nun 

105   mm 

155   mm 

A 

k  = 

A 

e   = 

1360.498 
0.871 

1991.151 
0.769 

2349.845 
0.359 

Ine  followinc  figure  displays  the  conditional  probability  of 
suppression  with  the  atove  estimators  for  the  parameters. 


P(S/9,r)  =  exp 


L  x 


r  (1  -  e  cos  9) 


] 


for    8=0    i.e.    along    the    main    direction    of    sight. 
P(S/e      =  0°,r)    =   exp     -  -  r    (1  -   e) 


(3.40) 


(3.41) 
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Range  of     r  i 
P  : 


•100     200 
0  1 


81  mm 

K  -  1360.498 
Ep3  -         0.871 


:  p 


105  rnm 

K  =  1991.151 
Eps  =■        0.769 


•  P 


155  mm 

K  -  2349.  845 

Eps  -        0.859 


Figure   15  -  FUNCTION        P(S/8  =  0°,r)   =•  cxp[-  £  r2(l   -  c)] 
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For  the  above  selected  values  of  the  standard  deviation  a, 


the      fire      rate 


f  ' 


and      the    parameters    K,     e    ,    and   H   the 


values   of    the    probabilities 

K 


P(s)  = 


(3.41) 


2a 


1  + 


-  e 


2a' 


PK(S)=  P(S)   - 


K 


-    G 


(3.42) 


POO  = 


a   /H  +  1 


(3.43) 


and   rates 


X    (wiLh  no   killing)    =  \      P(S) 


X    (with  killing)    =  Xr   P„(S) 
s  f      K 


(3.44) 


X   =  X   P(K) 
K 

are   displayed   on   the  succeeding  two  figures  cumber  16  and 

17.   Numerical  values  for  the   rate   of   rise   ^    are   also 

u 

computed   and   displayed  in   figure   17.  The  collected  data 

(suppression   intervals)  were    not    distinguished    among 

individuals   as   in   the  course   of  tha  foregoing  analysis, 

aecause  of  this  fact  the  numerical  values  fcr   X   are   based 

u 

en   the    expectation   of   duration   time   with   regard   to 
population. 


E[t]  = 


XfP(S) 


X  P(S)t 
e        fpCOdt 


(3.45) 


It   might   be   worthwhile   to   evaluate  the  rate  of  rise 

by  applying  the  Gamma   distributions   for   f  (t)   which   were 

discussed   in   chapter   III.E.2,   and   contrast   them  to  the 

values  of   X   gained  by  the  data  itself. 
u  3       1 


89 


Hence  solving  equation  (3.45)  by  having  set 


fp(0  =  ?^  (^)r_1  x 


(3.46) 


we  find: 


E[t]  = 


XfP(S) 


-  XfP(S)t  -Xt 
/   e        It^t  (Xt)r-i  X  dt  -  I 
0 


r(r) 


(3.47) 


if       X>  X  P(S)    and    by    changing   the   variaole    cf   integration: 


(X  -  XfP(S))t  =  z 


(X  -  XfP(S))dt  =  dz 


(3.48) 


we   get: 
EM 


-z     r-1 

e       z 


XfP(S)[X    -  XfP(S)T    0 


_      r    e z 

r  J      r(r 


dz   - 


XfP(S) 


=  1 


E[t]    = 


XfP(S)[X  -  xfp(s)]r      Afp(s) 


(3.49) 


Ey  using  the  earlier  developed  estimators,  we  receive  the 
mathematical  expression  for  the  rate  cf  rise  based  on  a 
Gamma  distribution  for  the  suppression  interval  data. 

r  r 


X  = 
u 


LXfP(S)[X  -  xfp(s)]r 


^fP(S) 


-1 


(3.50) 


3emarJc:  Ihis  derivation  is  only  true  for  A  >  X  P(S) ,  which  means 
that  for  applying  this  formula,  the  fire  rate  X-  may  not 
exceed  the  value: 


X,  < 


(3.51) 


f  -  PCs) 

Otherwise   the  integral    (3.47)     explodes  towards  infinite 
and  the  resulting  rate  of  rise  would  oe 


X  =  0 
u 


(3.52) 
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which  is  consistent  within  the  set  of  assumptions. 

Instead  of  using  the  Gamma  distribution  the  data 
themselves  were  applied  in  a  second  step  to  calculate  an 
estimate  foi  E[t]  and  ^u  based  on  the  number  of 
observations  n. 


E[t]  = 


XfP(S) 


\n  1=1 


n    XfP(S)t±\ 


-  1 


,  (3.53) 


Ihe  approximate  value  for  the  rate  of  rise  is  then 


X  =  - 

U   E[t] 


(3.54) 


If  killing  as  an  additional  event  is  considered,  the  rate  of 

rise  is  computed  with  the  same  formulas   (3.49)   and   (3.53) 

except   that   this   time   P..(S)  is  used  instead  of  P(S)  .   The 

K 

computation  of  the  values  for  the  formulas  (3.50  )  and  {  3.54) 
was  performed  by  APL  and  FORTRAN  respectively.  The  programs 
are  displayed  on  appendix  K. 
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PROBABILITIES      P(S),    P    (S) ,    AND     P(K) 


81  mm 

P      ^"\. 

25   m 

30   m 

35   m 

P(S) 

0.573 

0.495 

0.4  30 

Pk(S) 

0.503 

0.446 

0.394 

P(K) 

0.137 

0.100 

0.075 

105   mm 

P                ^\^ 

20  m 

25   m 

30   m 

p(s) 

0.731 

0.643 

0.564 

Pk(S) 

0.582 

0.542 

0.492 

P(K) 

0.272 

0.193 

0.142 

155  m 
P              N. 

30  m 

35   m 

40   m 

PCS) 

0.573 

0.497 

0.432 

Pk(S) 

0.480 

0.428 

0.378 

P(K) 

0.181 

0.140 

0.111 

Figure     16    -       NUMERICAL    VALUES    FORP(S),    P    (S)  ,    AND     P(K) 
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Row  (4)  and  (6)  and  row  (5)  and  (7)  of  figure  17  display 
comparable  values  for  the  rate  of  rise.  As  it  was  stated 
earlier,  the  Gamma  distribution  underestimates  the  related 
frequency  of  the  data  and  hence  the  comparable  values  for 
the  rate  of  rise  Au  in  row  (4)  are  larger  than  the  values  in 
row  (6).  The  same  is  true  for  row  (5)  and  (7).  They  differ 
in  general  by  10  -  M%.  If  we  are  willing  to  live  with  this 
fact,  the  rates  of  rise  evaluated  by  the  fitted  Gamma 
distribution  are  a  good  approximation  for  the  values 
computed  by  the  data. 

In  order  to  verify  the  four  Lanchester  equations  (2.76), 
(2.77),  (2.96),  and  (2.99),  presented  in  chapter  II. C,  it 
is  necessary  to  specify  particular  rates  given  in  the  figure 
before.  In  case  of  a  certain  known  composition  of  fire 
units,  specific  rates  could  be  developed  as  inputs  fcr  the 
model  equations. 

By  this,  the  model  equations  receive  their  specific  shape 
and  scale,  their  general  behavior  remains  the  same,  as  can 
be  seen  in  figure  9  and  10. 
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IV.    DISCUSSION 

Eased  on  the  research  effort  of  this  paper  it  was  found 
that  the  phenomenon  of  suppression  as  defined  earlier',  is  a 
multidimensional  problem.  It  is  influenced  by 
psychological,  physiological,  and  environmental  variables. 
Ihere  exist  many  possibilities  to  model  dependencies  in 
general  form,  however  to  make  quantitative  statements  about 
suppression  the  modellers  have  to  restrict  their  efforts  to 
those  variables  which  are  observable  and  measurable.  Since 
the  main  objective  of  this  paper  has  been  to  establish 
models  which  are  able  to  express  relationships 
quantitatively,  the  main  thrust  has  bean  to  formulate 
suppression  as  functions  of  weapon  systems  and  their 
dispersions. 

For  the  evaluations  of  the  models,  a  set  of  simplifying 
assumptions  was  necessary  in  order  to  guarantee  a 
mathematical  transparency.  Ihe  dependencies  developed  in 
this  thesis  postulate  some  satisfying  results  in  modelling 
suppression.  Ihe  models  reflect  sufficient  accuracy  of: 

The  physionomy  of  the  human  being  and  its  resulting 
behavior  with  regard  to  suppression.  Suppression  is 
mainly  caused  by  visual  and  auditory  perception. 

Ihe  influence  of  the  weapon  systems  i.e.  their  size 
and  their  firing  capability  are  determining  factors 
for  the  amount  of  suppression. 

Cf  course  the  detailed  results,  i.e.  the  estimation  of  model 
parameters,  are  based  on  selected  weapon  systems  and 
scenarios.  Many  possibilities  exist  for  further  work  in 
this  area  particularly  under  the  aspect  of  including 
extended   human   factors   components   in  the  construction  of 
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suppression  models. 

The   authors   feel   the  paper  may  provide  a  contribution  to 

future  design  of  vargames  and  simulations  as  well  as   weapon 

systems.    It   also   supports  a  more  careful  analysis  cf  the 
combat  situation. 
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APPENDIX  A 


ORIGINAL  DATA  81  mm 
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APPENDIX  A 


ORIGINAL  DATA   105  mm 
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APPENDIX  A 

ORIGINAL  DATA    155am 


f-i  f>    r-»  CS    h    rl    W    H  C<  «-«  «-* 


(N     CNC*<N(N<NCNCNCN<N<NCN(NCNCN(NC<<MCNCN* 


O    OOr-tOOOOOOOOOOOOOOOO 


ooooooor,<N(0-j-cN(Nm(jtcNuia>j3- 
r»  r-i0r-^or'^-.or-'r-f0ocN00<NO^«T-4 

«-«    (N  r(    CN  h    (N  r-tCN  .H    CN  H    OJ    rl    H 


o   oooooco<Nr*03<Nr*nnj   «  n  J  n  rt 


«-*  c*  o   ^vntor*-co<7»o^«<Nn3"i/»or*cocT»o 


99 


APPENDIX  A 


SUPPRESSION   IWTiKVALS  IN  SiiCONDS 
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APPENDIX  A 


105  m  SUPPttSSSION   INTERVALS  IN  SECONDS 
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APPENDIX  B 


DERIVATION  OF  CONSTANT  A 


The  constant  a  can  be  determined  corresp ondingly  to  the 
following  derivation: 

Modifying  equation  (3.  10)  in  Cox,   The   Analysis   Of   Einary 
Data,  p.  33,  we  can  define  a  transform  as: 


Zn 


S.  +  a' 


(1) 


Starting  with  the  original  model 

f(6,,r.) 


P.  =  Pi(S/9,r)  =  e 


l  l 


and   using   the   log- transformation 

y .   =  Zn  p . 
i  i 

wa   subtract 

/3i+a\ 

yi  -  yi =  in  \~nT~y "  "n  Pi 

and  choose  the  constant  a  such  that 

z[y±  -  y±]  =  o 

Since  s-Einoaial   (n  p   ,  "iP-^1  ~  pi^ 
which   can    be    approximated    by 


S     =  n.P(S/9,r)    +  /n.    U 
i         ^  * 

where   U    is   a   random   variable    with 

E[U]   =  0 


E[IT]  -  P1(l-Pi) 
Inserting  S  in  equation  (3) 


(2) 


(3) 
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Ely,  -  yt]    -  E 


in 


S^  +  a 


-  *n  P, 
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n.p .    +  /n7  U  +  a 
1    i i 
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An      1  + 
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/ — '  n.p. 

i   ri 


Approximating   this   by    Taylor    series 


=  E 
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i   Pi 
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31  /!-•  _  n.p 


n.p.  11/ 

-     11  ' 

It   is   sufficient   to   consider   the  firs-  two  terms  of  the 

Taylor  expansion  since  U  and  a  are  small  in  comparison  to  n. 


=  E 


JL_+      a 
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Ua 


/ — '  n.p.        2n .  p .  3/22        „    2  2 

'n.p.  ii  li        n .      p .         2n . p . 


11  ii  ii-1 

This    equation    has    to    be    zero    according    to     (3) 

Hence 


terms  cf  crder 


E[U]  =  0 


terms  of  crd  =  r    — 

n. 

i 


i   2p 


i  J 


a  -  -j  (1-  PjL)  =  0 


a  «f  <1-Pl) 


terms  of  higher  order  are  neglected. 
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APPENDIX  C 


DEBIVATION  OF  VARIANCE 


The  weighting  factor  Var  (In  f)  for  method  7  in  chapter 
III  B  is  computed  as  fellows: 

Expressing  the  conditional  probability  of  suppression 
E  (S/9,r)  as  a  guotient  of  the  nuinDer  of  successes  Si  and  the 
number  of  trials  nt  leads  to: 

i   n. 

1 


where  the   expectation   cf   f    is: 
Eff.J   =  P.(S/8,r) 


in  f .   =  In  \  — 
Vn. 


Var[£n   f .]    =   Var 


In 


l  -1 


Since   the    random   number    S   is    binomial    with    mean 
E(S)    =  nP(S/Q,r) 

and   variance 

Var(S)   =  nP(S/0,r)-(l  -  P(S/0,r)) 

it  can  be  stated  as  a  function  of  a  random  variable  U      with 
mean 

E(U)  =  0 

and  variance 
Var(U)  =  1 
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i.e. 


}   =  n.-Pi(S/0,r)  +  /i^P^S/G  ,  r)  •  [1  -  P.(S/G,r)]-U 


Inserting  this  in  the  equation  above: 


Var[£n   f.]    =  Var 


Var[£n   f . ]    =   Var 

x 


Var  [Jin  f .  ]   =  Var 


In 


n.P    (S/0,r)   +  /n.P. (S/0,r) [1   -   P . (S/0 ,r) ].  U  ' 


In      P. (s/6 ,r)-  I   1   + 


'[1   -   P.(S/0,r)] 

n.P. (S/0, r) 

x  x 


[1   -   P.(S/8,r)] 
J,.   P.(S/9,r)    +  JLn(l  +         n   P    (S/9>r)-   •    i; 


i   l 


Since   P(S/9,r)    is    constant    we   know   that 

VarUn   P(S/0,r))    =   0 
which    leads    to: 


Var[2n   f . ]    =   Var 


la      1  + 


'1   -   Pi(S/0,r) 

n.P.  (S/0, r) 
i   i 


knowing    that    for   small   x1 s 

£n(l  +  x)   ~  x        for      x  << 
we    can    state 


Var  [Jin   f ']  ss  Var 


1  -   P.(S/0,r) 

n.    P.(S/0,r) 
l      l 
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APPENDIX    D 


ANALYSIS    OF    RESIDUALS    81   mm 
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21*9. 1*55 3i 

EPS 
9.8538*33369 

SCAT   A 
RABGF   OF    It    ~i       0 
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0  .220*88229* 

"l  ,  77301033  7 
0.  J977f.  J7W3 

"C.  36517073  f,  7  6 
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ANAUSIS  OF  TOE  RESIDUALS  155  ma 
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APPENDIX    G 


ANALYSIS    OF    RESIDUALS    (METHODS) 
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APPENDIX    H 


ANALYSIS    OF    PARAMETERS    K     and      e 


Parameter  K 


DATA    81//// 


DATA   105MM 


DATA    15  5///-/ 


Method  5 


SCAT   KY 
A'l'S  OF  Xl    0       8 

a-::;  or  t:   1577.9     1573.3 


SCAT    KY 
■Anr.r  or  Xt    a     t 
i.lr.r.    OF   J\    J0K.7S       JOK.SS 


•  ♦  •         • 
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SCAT    KY 
.V.'Tf    *r   X:    0       8 
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5C^7    XY 

?Ar:cn  or  Xt  o     a 
ta/tg/:  or  it  2300     2550 


«     «       • 
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ANALYSIS   OF  PARAMETERS      K     AND      e 


Parameter  £ 
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ANALYSIS   OF   PARAMETERS      K     AND      e 
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APPENDIX    I 


API    REGRESSION    PROGRAMS 


NON-ITERATION  METHODS 

V  BE  ORES 

Cl]  Z«-    36    2    pi 

T2]  Zlill+U 

[3]  Z[;2]*-7 

[4]  zt-<-§z 

[5-]  BETA-{  (HCZT+.xZ)  )  +  .*Z?)+.*Y 

[6]  K*Cl)*SETAZll 

[7]  £PS«-#xfl£T/tn2  3 

[  8  ]  YHA  T+  (  n.?T.4  [  1  3 x  £/ )  +  (  Bff2M  [23*7) 

[9]  RES+Y-YHAT 

[10]  /I*    36    2    pi 

ni]  d C  :  1 3+J 

[12]  At  ;2]*-RES 


ITERATION  METHODS 

V  ITER    P 

[I]  ST<7-«-l 

[2]  7«-(©(  (/?+(0.5*(l-P)  )  )* ,v)  )v5I<7 

[3]  U+UOiSIG 

[4]  y+-70vSIG 

[5]  Z-*-    36     2    pi 

[6]  Z[;l]«-0 

[7]  ZT;2>K 

[8]  ZT*-§Z 

[9]  BFTIT^i  (SKZT+.xZ))  +  .xZ2,)  +  .xy 

[io]  k«-(~i  hbhtitTi  3 

[II]  EP5<-K*BET 'It T2  D 

[12]  YRAT+(BETIT[1  ]*UO)  +  ( BETITl2~\*VO) 

[13  3  RES+-YC-YHAT 

[143  /I*    36    2    pi 

[15  3  4[;13-YC 

[16  3  4[;2>-ffff5 


REMARK:  GIG  CHANGES  IN  ACCORANCE  TO  THE 
DIFFERENT  ITERATIVE  METHODS. 
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APPENDIX  J 
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APPENDIX  K 


PROGRAMS  FOR  COMPUTING  RATES  OF  RISE 


Program  for  formula  (3.54) 


FORTRAN  IV  G  LEVEL   21 


MAIN 


0001 

DIMENSION  TIM(500)  rSLA(20) ,  RETAU(20) 

0002 

N=348 

0003 

READ  (5,10)   (SLA(I) ,1=1,18) 

0004 

10 

FORMAT   (F10.5) 

0005 

READ(5,11) (TIM(I) ,I=1,N) 

0006 

11 

FORMAT (F10. 5) 

0007 

DO  1   1=1,18 

0008 

SUM=0 . 0 

0009 

DO  2   J=1,N 

0010 

C=EXP( (SLA(I) *TIM(J) ) ) 

0011 

SUM=SUM+C 

0012 

2 

CONTINUE 

0013 

A=SUM/FLOAT(N) 

0014 

B=A-1.0 

0015 

ETAU=B/SLA ( I ) 

0016 

RETAU(I)=1.0/ETAU 

0017 

1 

CONTINUE 

0018 

WRITE(6,12) (RETAU(I) ,1=1.18) 

0019 

12 

FORMAT(1X,3F10.5) 

0020 

STOP 

0021 

END 

Program  for   formula    (3.50) 


V     LAMV 
[  1  ]  LA+  (LAMS*    ( {LAM-LAMS  >#  )  H  (LAM*R  ) 

[2]  L4£AKU((1tL4  )-(ULAMS)) 
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